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DMITRII VASIL'EVICH EFREMOV 


Professor Dmitrii Vasil’evich Efremov, member of the CPSU, died on November 27, 1960 in his 61st year, after 
a long and serious illness. 


D. V. Efremov was born in Leningrad in 1900. After graduating from the Leningrad Polytechnical Institute, 
D. V. Efremov worked at the ‘Elektrosila' factory and rose from an ordinary engineer to become the chief engineer. 


During the last 14 years,D, V. Efremov was the Deputy Minister and then Minister of Electrical Industry of the 
USSR, First Deputy Minister of Electrical Power Stations and Electrical Industry, Deputy Chairman of the Bureau for 
Electrical Power, Chemical and Forestry Industry of the Council of Ministers of the USSR, and Deputy Director of the 
Atomic Energy Utilization Directorate of the Council of Ministers of the USSR. 


D, V. Efremov directed the construction of major Soviet turbogenerators and many electrical machines, as well 
as special nuclear particle accelerators. 


D. V. Efremov, a scientist and communist, was a great patriot and continued his work right up to the last days 
of his life. 


During the 19th congress of the Party,he was elected as a candidate for the Central Committee CPSU. 


For his distinguished services to his country, D. V. Efremov was awarded three orders of Lenin, three Labor Red 
Banner orders, a ‘Mark of Distinction’ order, and other medals. He was awarded the Lenin and Stalin prizes. 


D. V. Efremov was a pleasant and considerate person and was very concerned with the welfare of his colleagues. 


The memory of Dmitrii Vasil'evich Efremov will always remain in our hearts. 


Pravda, November 29, 1960 A group of colleagues 


989 


235 
NEUTRON SPECTRA EMITTED IN THE FISSION OF U 


AT 0, 45, AND 90° TO THE DIRECTION OF MOTION 


OF THE FRAGMENTS 


Yu. A. Vasilev, Yu. S. Zamyatnin, E. I. Sirotinin, 
and E. F. Fomushkin 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 449-4054, 
December, 1960 
Original article submitted February 29, 1960 


Neutron spectra emitted at 0, 45,and 90° to the direction of motion of the fragments in the 
fission of U* by 14.3 Mev neutrons have been measured. The angular distribution of the 
fission neutrons was found to be N(0°) : N(45°) ; N(90°) = (3.23 0.12): (1.75 + 0.07): 1.00. 
The experimental data are compared with spectra calculated on the assumption of an iso- 
tropic Maxwellianneutron distribution in the center-of-mass system of the fragment. 


Introduction 

Measurements of fission neutron spectra [1, 2] and the angular distribution relative to the direction of 
motion of the fragments [3, 4] have shown that the main mechanism responsible for the production of fission neu- 
trons is evaporation from the excited fragments. Spectra calculated on the assumption that the neutrons are evap- 
orated isotropically and have a Maxwellian distribution in the center-of-mass system of the fragment [1, 5] are in 
good agreement with experiment,in spite of the simplifications introduced into the theory. However, the values 
of W (the average kinetic energy per nucleon in the fragments) which are in the best agreement with the experi- 
mental data on neutron spectra are always lower than the directly measured values [6, 7]. This leads to some 
doubt as to the correctness of the assumptions made in the theory. 


Complete information on the mechanism of production of fission neutrons can be obtained from studies 
of the spectrum and angular distribution of the neutrons in the center-of-mass system of the fragment for various 
cases of fission. However, it is relatively difficult to obtain such data from the spectra of fission neutrons ob- 
tained in the laboratory system of coordinates,since it is necessary to allow for the velocity distribution of the 
fragments, their excitation energy distribution [8], and finally, the averaging of the spectra over all angles rela- 
tive to the direction of motion of the fragments. It is very difficult to set up an experiment in which all these 
effects would be excluded, and hence in the first approximation, it is expedient to measure neutron spectra emit- 
ted at some definite angles 6 to the direction of motion of the fragments. This excludes the necessity of averag- 
ing the spectra over the various angles, and thus considerably simplifies the transition from the laboratory to the 
center-of -mass system. 


In the present work,we have determined the spectra of neutrons emitted in the fission of U® by 14.3 Mev 
neutrons at 0, 45,and 90° relative to the direction of motion of the fragments. The angular distribution of the 
neutrons was thus obtained at the same time. 


Measurements and Results 
The spectra and the angular distribution of the neutrons were measured by the time -of- flight method. 


The flight-path was 75 cm and the resolving time 7-107*sec. A detailed description of the method and appara- 
tus is given in [9]. 


The particular feature of the experiment was that the source of the fission neutrons was a multilayer fis- 
sion chamber with fragment collimation (Fig. 1). 


The direction of motion of the fragments in this chamber was defined by collimators made from nickel 
ribbons 0.1 mm thick. Each collimator cell was in the form of a rectangular parallelepiped with a base of 
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2x 2mm and a height of 4mm, The average and most probable fragment exit angles were roughly 14°. This col- 
limation enabled us to record up to 1.2% of the fissile events. 


The U* films were deposited on either side of aluminum discs 0.5 mm thick, The diameter of the films was 
80 mm and their thickness about 6 mg/cm”. The total weight of the U™ was 3.5 g,and 0.75% of the fissile events 


occurring in the uranium were recorded during the experiments. 


The chamber was filled with a mixture of argon and carbon dioxide (10%) to a pressure of 760 mm Hg. With 
the high-voltage electrodes at about 1000 volts, the rise-time of pulses from the chamber was approximately 0.1psec. 


The position of the chamber relative to the target 
of the accelerator tube is shown in Fig. 2, The measure- 
ments were carried out at 0, 45,and 90° in the neutron 
energy interval 0.4-5 Mev. About 2.5 x 10° fission neu- 
trons were recorded at each angle. The background was 
of the order of 10 to 15% and increased with increasing 
angle. At the limits of the interval, i.e., at 9 = 90° the 
background reached 45%, Statistical errors did not ex- 
ceed 10-15%, except for the energy range 0.4-0.6 Mev, 
where they reached 25% at 90°. 


Figure 3 shows the neutron spectra emitted in the 
fission of U™ at 0, 45,and 90°, Figure 4 shows a compari - 
son between the estimated integral spectrum obtained on 
the basis of the present results with the measured spectrum 
given in [9]. The°two sets of results agree to within ex- 
perimental error. 


Fig. 1. Fission chamber with the fragment collimation 


system. 

The neutron spectra emitted by the fragments were 
determined from the experimental data by subtracting the 
neutron spectra evaporated prior to fission. This was done 

Deuteron beam with the aid of the results given in[9]. It was assumed 
Target that the angular distribution of the neutrons evaporated 
- prior to fission is isotropic. The spectra obtained after 
g this subtraction are shown in Fig. 5. 
oa PeAtattOntal ——— SS 5 The angular distribution of the fission neutrons 
U eect Se aaa 5 (cf. Table) was determined by numerical integration of 
deposits aHHSHHLHARLE WH 90° the spectra shown in Figs. 3 and 5 with respect to the 
a ee 
pn 7 (ate Dissent data wets also used oe SAanENe the 
dua MVUUERAEHLACAUGERLUDEL angular distribution of fission y -rays with enere\cs Ey? 
TOO TIONOCUOGBDTEAOL > 0.3 Mev. Assuming that the y spectrum is indepen- 


dent of the angle of emission of the y rays relative to 
45? the direction of motion of the fragments, it was found 
that ny (0°) ; ny (45°) : ny (90°) = (1.31 + 0.07) ; (1.22 4 
+ 0,06): 1.00. 


ie 


Fig. 2. Position of the fission chamber relative to the 

target of the accelerator tube. The arrows show the Diecussion of Results 

directions to the neutron detector. As was to be expected,the present results indicate 

considerable differences between neutron spectra (in the 

laboratory system) emitted at different angles 9 to the direction of motion of the fragments. There is a considerable 
anisotropy in the angular distribution of the neutrons : by= N(0°)/ N(90°) = 3.23 + 0.12. The result b'y = 4.03 + 0.23, 
which was obtained after subtraction of the neutrons evaporated prior to fission, agrees to within experimental error 
with the corresponding quantity bT = 4.35 + 0.19 obtained for the fission of Wie by thermal neutrons [3]. The some- 
what lower value obtained for b{4 as compared with bt may be explained by an increase in the average energy of 
neutrons evaporated from the fragments [9]. The agreement between the values obtained for b‘y and bT confirms the 
value obtained in [9] for the fraction of neutrons evaporated prior to fission. 
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Since the form of the neutron spectra at different angles 9 and the magnitude of the anisotropy are pee depen- 
dent on the neutron spectrum in the centerof-mass system of the fragment, we have tried to deduce fede Es on 
this spectrum from the results of our experiments. With this aim in view, we have carried out calculations using va- 
rious assumptions about the neutron spectrum in the center-of-mass system of the fragment. In these calculations we 
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Fig, 3. Neutron spectra in the fission of ee by 14.3 Mev neutrons at 0° ®. 


45° P), and 90° ch) relative to the direction of motion of the fragments. The 


errors are statistical and all spectra are reduced to the same number of fissile events. 
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Fig. 4. Comparison of the estimated integral spectra based on the present results ©) 


with the spectrum reported in [9] (continuous curve and 9). The errors are statistical. 
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took into account the fragment mass distribution, the dependence of the total kinetic energy of the fragments on Rg, 
i.e., the fragment mass ratio [7], the dependence of v on the fragment mass [10], and also the angular distribution of 


fragments passing through the collimator. In all cases an isotropic neutron emission in the center-of-mass system of 
the fragment was assumed, 


Angular Distribution of the Fission Neutrons of Uae 


Measured Calculated 
neutrons from spect 
all neutrons {fragments T=1 Mev fF =0.6 Mev | T=0,2 Mev on [12] 


N (0°) 3,2340,12 | 4,03-0,23 1,75 3,82 42,40 4,30 
N (45°) 4,75--0,07 1,890, 12 4,54 2,88 17,10 2,42 
N (90°) 1,00 1,00 4,00 1,00 1,00 1,00 


Note. The calculated angular distributions of the neutrons emitted by fragments were obtained on the 
assumption of a Maxwellian neutron distribution in the center mass system of the fragment, 
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Fig. 5, Neutron spectra emitted by the fission fragments of UT at 0° (O), 45° (@), 
g Pp y g 


F(E), arb, units 


and 90° ch) to the direction of motion of the fragments. The errors are statistical. 


The results of these calculations for a Maxwellian spectrum in the center-of-mass system, corresponding to 
temperatures of 0.2, 0.6,and 1,0 Mev, and the spectrum obtained in [11] for Cf’, are compared with experimental 
data in Fig. 6. As can be seen, the experimental results agree neither with the calculations nor with the spectrum 
given in [11]. 


Calculations are being carried out at the present time which will take into account the temperature distribu- 
tion as described in [12]. According to our estimates,this should increase the fractionof low- and high-energy neut- 
rons in the spectrum, and should, at any rate,considerably reduce the discrepancy between theory and experiment. 
The discrepancy between the calculated and experimental angular distributions (cf. Table) may possibly be due toan 
anisotropic emission of neutrons in the center-of-mass system of the fragment [13], If it is assumed that the preferred 
emission of neutrons in the backward and forward directions is symmetric with respect to the angle of 90° [12, 14],then 
corrections for the anisotropy would change the form of the spectrum in the same direction as the temperature distri- 


bution. 
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It should be noted that the apparent discontinuity at about 0.6 Mev for 9 = eOgsneule not be crestes a se 
ously, because of the inadequate accuracy in this region. However, if this peak soese fact occur, then it mig 
explained by the emission of a small number of neutrons from the waist of the dividing nucleus. 


500 


400 


300 


F(E), arb. units 


200 


100 


300 


200 


F(E), arb. units 


100 


200 


100 


F(E), arb. units 


E , Mev 


Fig. 6. Comparison of experimental results () with neutron spectra calculated on the assumption of 


4 Maxwellian distribution of neutrons in the center-of-mass system of the fragment with T = 0.2 Mev 


(————), 0.6 Mev (— - —), 1.0 Mev (— - - —) with the spectrum obtained in [11] (— — —), 
All the spectra are normalized at @ = 0°, 


In conclusion the authors wish to thank P. V. Toropov for assistance in the preparation of the electronics; Yu. 
Ya. Glazunov, A. N. Maslov, N. I, Nemudrov, V. A. Parshin,and V. S. Khorkhordin for taking part in the measure- 
ments; and V. A. Komarov, M, P. Novikov, G. A. Peretokin,and L. A. Chernov for writing the computer program. 


994 


LITERATURE CITED 


B, Watt, Phys. Rev. 87, 1037 (1952). 

B. G. Erozolimskii, Atomnaya Energiya, Supplement No. 1, p. 74 (Atomizdat, Moscow 1957).° 

J. Fraser, Phys. Rev. 88, 536 (1952), 

R, Wilson, Phys. Rev. 72, 189 (1947); R. Ramanna and P, Rama Rao, Communication No, 1663 submitted by 
India to the 2nd International Conference on the Peaceful Uses of Atomic Energy (Geneva, 1958). 

5. B. G. Erozolimskii, Atomnaya Energiya, Supplement No. 1, p. 74 (Atomizdat, Moscow 1957).* 

6. D. Brunton and G. Hanna, Canad. J. Research 28A, 190 (1952); R. Leachman, Phys. Rev. 87, 444 (1952). 
is W. Stein, Phys. Rev. 108, 94 (1957). 
8 
9 


Ao 


R, Leachman, Phys. Rev. 101, 1005 (1956); R. Leachman.and S. Kazek, Phys. Rev. 105, 1511 (1957). 
; Yu. A. Vasil'ev et al., Zh. @ksperim. i teor. fiz. 38, 671 (1960). 

10. S. Whetstone, Phys. Rev. 114, 581 (1959); V. F. Apalin et al., Atomnaya Energiya 8, No. 1, 15 (1960).° 

11. H. Bowman and S, Thompson, Communication No. 652 submitted by the USA to the 2nd International Con- 
ference on the Peaceful Uses of Atomic Energy (Geneva 1958). 

12. J. Terrell, Phys. Rev. 113, 527 (1959). 

13. OD. Hill and J. Wheeler, Usp. fiz. Nauk 52, Nos. 1 and 2 (1994). 

14. ‘V. M. Strutinskii, Nuclear Reactions at Low and Intermediate Energies, a collection of papers [in Russian] 
p. 522 (Academy of Sciences Press, Moscow, 1958). 


* Original Russian pagination, See C, B, translation, 


995 


COHERENT ELECTRON RADIATION IN A SYNCHROTRON 


Yu. M. Adovand | V.V. Elyan 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 455-460, 
December, 1960 
Original article submitted May 28, 1960 


The coherent radiation of electrons in a synchrotron was investigated experimentally for a wavelength 
of 10 cm, corresponding to the 50th harmonic of the electron frequency of revolution for various elec- 
tron distributions inside the separatrix. In most cases, the theory of coherent radiation of electrons is 
confirmed, A disparity with the theory is observed for almost uniform electron distributions. The radia- 
tion was used for measurements of the electron phase oscillations and adiabatic damping of the phase 
oscillation amplitudes, The experiment was carried out on the synchrotron at the Physics Institute, 
Academy of Sciences, USSR at 280 Mev with an initial betratron acceleration. 


Introduction 

In the case of circular motion of an electron, electromagnetic waves are radiated with frequencies that are 
multiples of the frequency of revolution [1]. In the region of wavelengths comparable with the dimensions of the 
bunch, the radiation is coherent and depends on the character of the electron distribution around the orbit[2-6]. This 
is taken into account by the introduction of a form factor Fnm in the formula for the coherent radiation intensity 


Wom = CO Loras (1) 


where Wnm is the intensity of the radiation of an individual electron undergoing phase oscillations with a frequency 
2; N is the number of accelerated electrons; n is the number of the harmonic; m = 0, 1, 2, 3, . . . characterizes the 
side frequencies close to the nth harmonic. 


The form factor Fnm has been calculated for several shapes of the stationary distribution of the electrons for 
the phase oscillation amplitudes f (#) [5): 


for 
ft (®) = x5, V 20,7 — 
20) 0 y (2) 
_ Sin? (2ny) 
no (2n@,)? ’ (25 
for 
f(D) =8(9—Q,), (3) 
die = i (n®,); (3') 
for 
my, d\2 
f(O) = 20 o7 (ae) 
ie 
(4) 
nr? 
Hee = rig oe ; 
(4") 
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where ®p is the mean amplitude of the phase oscillations in the bunch. 
If m ~# 0, then Fpm = 0 for all three cases of the stationary distribution. 


The acceleration of electrons in the synchrotron is accompanied by the damping of the phase oscillations. The 
coherent radiation intensity then changes in accordance with the laws described by the form factors (2°), (3°), and 
(4"), depending on the character of the electron distribution over ®. In cases (2") and (3'), this change is of an os- 
cillatory character. This property was first observed on the synchrotron at the Physics Institute, Academy of Sciences, 
USSR at 5 Mev in a study of the coherent radiation for the 16th and 24th harmonics and was used for the measure- 
ment of adiabatic pinching of the bunch [6]. 


For a nonstationary phase oscillation amplitude dis- 
tribution of the electrons,it is not necessary that Fy; = 0 
for m # 0 [5]. Moreover, the radiation signal turns out to 
be amplitude-modulated by the phase oscillations frequen- 
cy. In this case, the frequency of the electron phase os- 
cillations can be measured experimentally. 


Equipment 

The radiation was brought out of the synchrotron 
chamber through a lateral tube of 8 cm diam, A superhe- 
terodyne receiver was used to pick up the radiation. The 
layout of the equipment is shown schematically in Fig. 1. 
A born antenna 1 and mixer 2 were placed right next to 
the accelerator. A klystron (K-11) oscillator 9, interme- 
diate-frequency amplifier 3, and stabilized power sup- 
plies 4 and 8 were placed at a distance of about 3 m from 
the accelerator magnetic to avoid effects from the vary- 
ing magnetic field on the operation of the circuit. The intermediate-frequency amplifier was six-stage tuned amp- 
lifier with a 0.5 Mc passband and amplification factor of ~5 x 10°. The intermediate frequency was chosen as 25 
Mc. After detection, the signal was fed to an OK-17M oscillograph 5 equipped with a photographic attachment. The 
oscillograph sweep was triggered by a pulse synchornized with the magnetic field of the accelerator. The oscillog- 
raph sweep was calibrated by time marks. The equipment allowed the receiver to be retuned over a range of ~10 
Mc, which was limited by the region of klystron oscillation. The klystron signal frequency was changed by the vary- 
ing of the voltage on the reflector by means of two selsyns 6 and 7. 


Fig. 1. Block diagram of the equipment. 


Experimental Results 
Oscillographs were taken for various distributions of the electrons in the separatrix. 


In the synchotron with betatron starting, the distribution could be changed by shifting the time of switching on 
the high-frequency accelerating voltage with respect to the optimum time. The experimentally obtained dependence 
of the relative number of particles N trapped in synchrotron operation on the instance for switching on the accelerat- 
ing field is shown in Fig. 2. The optimum time for switching on the accelerating field corresponds to N =100. On 
segments I and III, the electrons fall inside the separatrix on phase trajectories close to each other. Therefore, the 
character of the electron phase oscillation amplitude distribution is close to a 6 function. For region If, a more 
even filling of the separatrix and a ® distribution of the electrons close to (2) are characteristic, The amount of trap- 
ped electrons on segment I depends on the accelerating field voltage. Hence when this amount is measured, the fre- 


quency should be chosen from the greatest value of N on segment I. 


Figure 3 shows typical oscillograms of the radiation obtained when the accelerating field was switched on at 
segment II, As seen from the oscillograms, there was no radiation during the betatron operation of the accelerator 
as long as the electron energy was ES 4 Mev; it appeared only after formation of the bunch, At E > 180 Mev,the 
radiation intensity decreases sharply and signal drops below the level of the equipment noise. The shape of the ra- 
diation pulse depends on the conditions of trapping in synchrotron operation, i.e., on the amplitude and frequency of 
the accelerating field, its rise time, and the time at which the field is switched on, Under some conditions, the ra- 
diation appears to be partially amplitude-modulated, In the oscillogram of Fig. 3b, the region with this type of mo- 


dulation is indicated by an arrow, 
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Figure 4 shows oscillograms of this region taken with a fast sweep for three values of the accelerating voltage 
amplitude U. The modulation frequency ratio measured from the oscillograms 
By & Fo: F3 = 0.85: 1: 1.15 


is in good agreement with the theoretically calculated ratio of the phase oscillation frequencies: 


Fy: Fy: Fs3=VU,sing, 1 VU, 510 yy: VUssing,s =0,84:1: 1,10, 


where ¢g is the equilibrium phase. 


This shows that such amplitude modulation is due 
to the phase oscillations of the particles distributed under 
nonstationary conditions in accordance with the ampli- 
tudes of these oscillations. 


The table lists the frequencies of the phase oscil- 
lations F, the maximum amplitudes of the phase oscilla- 
tions max, obtained from the known relation between 
F and max [7], and the maximum size of the separatrix 
?max determined from the oscillograms of Fig. 4. 


| “le I eel Pea 
_ eee oe 


: eee ts 


100 150 200 250 300 


In the calculations we used the values E = 80 Mev, 


magnetic field exponent 0.65, frequency of revolution of 
Fig. 2, Experimental dependence of the number of par- __ the electrons 58 Mc. 


ticles N trapped in synchrotron conditions of acceleration 
on the time at which the high-frequency accelerating 
voltage is switched on. The time is counted from the 
instant at which the magnetic field passes through zero. 


aoe t, psec 


From a comparison of &max and J max,it is seen 
that the contribution to the modulation comes from the 
electrons grouped close to the separatrix. 


The radiation is of an oscillatory character only 


Data for the Phase Oscillation Frequency Measurements under specific conditions of bunch formation in the tran- 
sitional operation, when the amplitude distribution of the 
U, volts | 8!n 9 F, ke a max | ?°max phase oscillations is close to a 6 function. For example, 
when the accelerating field is switched on at times cor- 
800 0,77 52445 9249 | 100-440 responding to seo e I in Fig. = = radiation depends 
4000 0,87 6346 | 145+44 | 120-12 on the energy in the way shown in Fig, 5. Shortly after 
1200 0,94 84+6 | 127442 | 130413 the beginning of synchrotron operation, however, the ra- 


diation becomes either quite irregular (noise) (see, for 

example, Fig. 3a) or partially regular (Fig. 3c). Depend- 

ing on the frequency of the accelerating field, the regu- 

larity can begin to be disturbed at different electron ener- 
gies from E = 8 Mev (Fig. 5a) to E = 60 Mev (Fig. 5b). This effect is obviously connected with the redistribution of 
the electrons in the separatrix during the acceleration process, 


Slow oscillations of the radiation intensity (see Fig. 5) are due to the adiabatic damping of the phase oscilla- 
tions; similar oscillations have been observed in [6]. 


Figure 6 shows curves of the phase oscillation amplitude damping calculated by the method of adiabatic inva- 
riance [7] and the experimental results obtained on the basis of the analysis of the oscillograms (see Fig. 5). 


While the experimental points are in very good agreement with the calculated curves for ®) iS 140°, they are 
not in agreement with the curves for 140°S ©) i S$ 180°. This means that the capture of electrons in synchrotron oper- 
ation on segment I of Fig. 2 takes place on phase trajectories with amplitudes less than 140°, 


We also measured the bandwidth of the frequencies of the 50th harmonic. It turned out that, for radiation with 
a signal of oscillatory form (see Fig. 5), the bandwidth of the harmonic frequencies does not exceed 0.5 Mc. This 
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Fig. 4. Oscillograms of the intensity of coherent ra- 
diation amplitude- modulated by the phase oscillation 
frequency. A) U = 300 volts; b) U = 1000 volts; c) 
U = 1200 volts. 
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Fig. @. Modulation of the radiation intensity by slow oscillations. 


confirms the theoretical conclusion that for distributions of type (3) the side frequencies of the radiation vanish, i.e., 
Fnm = 0 for m # 0. 


The noise character of the radiation (see Fig. 3) is accompanied by a broad frequency band (Fig. 7), i.e., in 
this case Fnm * 0 for m # 0. This has not been explained theoretically thus far. We believe that the phenomenon 
of broad frequency bands near the harmonics of the radiation of the electron bunch with a distribution of type OO) can 
be properly explained qualitatively by the unavoidable fluctuations in the phase oscillation amplitude distribution of 
the electrons, which are not taken into account by the theory. 


The sharp decrease in the radiation intensity at 
electron energies above 180 Mev is explained by the fact 
that up to this time the phase oscillation amplitude dis- 
tribution of the electrons is close to Gaussian (4). The 
value of Fy) (4') then becomes incomparably smaller than 
(2) and (3). The formation of a distribution of type (4) 
is confirmed by other experiments made at the synchro- 
tron of the Physics Institute, Academy of Sciences, USSR 
at 280 Mev [8]. 


CONCLUSIONS 
1. By means of coherent radiation of electrons 
one can readily record the beginning of the formation 
of the electron bunch in a synchrotron with betatron start- 
ing. 


Fig. 6. Adiabatic damping of phase oscillations. 1 
and 2) Theoretical curves; O and @) experimental 
results for ®o; = 140° and ®oj = 180°, respectively. 


2. For a nonstationary distribution of electrons in 
the separatrix, the radiation can be used for the meas- 
urement of the phase oscillation frequency. 


3. The use of coherent radiation from the 50th 
harmonic for the study of phase oscillation damping is 
possible only for phase oscillation amplitude distributions 
of the 6-function type. For other types of distributions, 
lower harmonics should be used. 


4, The theory of coherent radiation of electrons 
in the synchrotron is in agreement with experimental only 
in special cases. To generalize the theoretical conclu- 
sions, it is necessary to consider the effect of the fluctua- 
tions in the electron distribution on the coherentness. 


“5 4 f,Mc 


Fig. 7. Relative distribution of the radiation intensity 
in the band of the 50th harmonic for optimum switch- 
ing on of the accelerating field. 


In conclusion, the authors thank Prof. P. A. Cheren- 
kov, Prof. M.S. Rabinovich, Prof. A. M. Prkhorov, and 
L. V. logansen for valuable discussions. 
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HEAT EXCHANGE DURING THE FLOW OF MERCURY AND WATER 
INGA ITIGHTEY PACKED ROD PILE 


V.. I. Subbotin, P. A. Ushakov, B. N. Gabrianovich, 
and A. B. Zhukov 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 461-469, 
December, 1960 
Original article submitted August 28, 1959 


The results of an experimental investigation of heat transfer in mercury and water in tightly packed 
rods are presented. The problems of thermal simulation of the rods are considered, the results of 
measuring the temperature field of the heat-exchange surface are given, and the experimental ap- 
paratus, procedure, and processing of the data are described. For a given group of tightly packed rods 
a procedure is suggested for making an approximate thermal calculation. 


It is customary in calculating heat transfer to liquids with Pr > 0.6 in complex channels torecommend the equa- 
tions derived for circular tubes, for instance the equation of M. A. Mikheev [1]: 


Nip 0,021 Rens Pre? (Pr iPra\s a, (1) 
the equation of Colburn [2]: 


Nu = 0,023Re%8Pr'/s, (2) 


and others. The hydraulic diameter dh is taken accordingly as the reference dimension. These recommendations are 
endorsed by a number of authors [3-6]. Recently, however, the results of investigations have been published showing 

that the use of dh as the defining dimension does automatically imply allowance for the peculiarities of heat transfer 

in rod piles subjected to an arbitrary flow pattern [7, 8]. Weisman [8] showed on the basis of a generalization of sev- 
eral papers that the influence of the spacing in a triangular latticework of rods can be taken into account by introduc- 
ing into Eq. (2) the following variable coefficient in place of the factor 0.023; 


K =0,026x — 0,006, (3) 


where x = s/ 2R, is the relative rod spacing and is equal to 1.1-1.5; Re > 20 » 10°, 


It is unfeasible to extrapolate the relation (3) in the hope of obtaining reliable data for closely packed piles 
with x = 1, since the heat exchange in closely packed piles depends on the construction and thermal conductivity of 
the rods, even more strongly than in piles with x > 1, due to the acute asymmetry of the cells. 


The published data on heat transfer of piles in axial flow of a liquid metal (Pr « 1) is almost completely re- 
stricted to the papers [9, 10]. The results of these papers, evidently, can be used for approximate heat exchanger 
calculations. It should be noted that in general no data have been published on tightly packed piles. 


Problems in the Thermal Simulation of Tightly Packed Rods 

Consider a pile with a large number of heat-generating rods. In this case the influence of edge effects on heat 
exchange and the hydrodyanmics of the central pile cells cannot be taken into account. We shall adopt the follow - 
ing assumptions: 1) The physical parameters of the system are constant; 2) the liquid is incompressible; 3) heat is li- 
berated uniformly throughout the volume of the rods; 4) ideal thermal contact is maintained between the heat-gene- 
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part of the rod and the protective casing; 5) the longitudinal cross convection of heat in the rods is small, 6) the tem- 


perature and velocity distributions are steady-state. 
As a result of analysis of the differential energy equations, the equation of motion, and appropriate boundary 


conditions on the basis of similarity theory a syste of dimensionless parameters is obtainable, characterizing the 


heat-exchange process in the pile cell: 
os l oat he P ) ; 
Nu=t(q Si? oe ee Re, fe Nae 
= l Ay Ag ) : 
sat (igi bd eerie 
(ie Saree Gs Q, ay ie ’ Ve , R 


aime ie ny = ab a f free (5) 


(4) 


where J is the pile length (m); € = R,/R, is the dimensionless radius (Ry, Ry are the inside and outside radii of the tube 
or rod casing, respectively, in m); 4, Xz, Ag are the heat conductivities of the heat-generating rod, rod casing (or 
tube), and liquid, respectively (kcal/m- Tit eC) aces fu he dX, is a dimensionless temperature parameter (tw is 
qk 

2 : i<} 
the wall temperature, °C; q the heat flow density, kcal/m + hr); t¢ is the temperature of the heat carrier (°C); dh = 
= 4F/P is the equivalent hydraulic diameter of the pile cell, m (P is wetted perimeter of the pile cell, m; F the ac- 
tive cross section of the cell). 


In the channels of the pile there is no axial symmetry, so that the temperature of the heat carrier and heat-ex- 
change surface is governed not only by the properties and flow conditions of the heat carrier, but also by the para- 
meters of the heat-generating rods ( £4, Xy/A2, A2/d3). 


In many instances it proves difficult to realize exact simulation of the rod, so that one must resort to approxi- 
mate simulation. Consider the two limiting cases: X; > Az and Ay *KAz. The conditions t = const or q = const will 
apply to the inner surface of the rod in correspondence with these two cases. 


Let us examine the conditions for approximate simulation in the case when A; Az. The heat-exchange pro- 
cess in the cell is uniquely determined by the geometry, flow conditions, properties of the heat carrier, and the boun- 
dary conditions, which in dimensionless form may be written 


i(zp.% Re, Br ie 11) <0, 
h q 


The connection between the parameters Ty and 1—? 


can be found by solving the problem of heat conduc- 


: . . q 
tion for the heat-generating rod and casing. If we write Tw in Fourier series form 


Ty = Dy Agp cos 6kq, (6) 
| 


the distribution of heat flow will have the form 


q—-4q es : = 
Epes Se By, cos 6kp = — SS 6k C6, Agr COS 6k; (7) 
=1 k=1 
se (on 
Pep adi EN erent ce iy (8) 


ce Met Ete m= es Wb . . 
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Fig. 1. Experimental heat exchanger: 1) Case; 2) flanges; 3) central tube; 4) 
side tubes; 5) rod displacers; 6) contour displacer; 7) thermocouple case; 8) pack- 
ing nuts for thermocouples to measure temperature of the mercury in the cells; 


9) thermocouple panel. 


In the first approximation in the expansions (6) and (7),we stop with the fundamental harmonic (k = 1). Then we can 


write 


Nu =i Re, Pr, ook 


(9) 
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re alte, Pr, a (10) 


The resultant system of similarity criteria for approximate simulation was used by the authors to generalize the 


experimental results. 


Experimental Apparatus 
The transfer of heat to mercury was studied on a stainless steel (type 1Kh18N9T) loop. The mercury was cir- 


culated by a centrifugal pump. The mass flow of mercury was regulated by a throttle device and was measured pre- 
liminarily by a graduated Venturi tube. 


The construction of the experimental heat exchanger is shown in Fig. 1. The dimensions of one heat-exchanger 
assemblies are the following (in millimeters); Ry = 8.879; Ry = 6.8 +0.025, 7 - 5005, 2, = 430**% dp = 1.81. 


The heat exchanger consists of a cylindrical case, inside which a pile of seven tubes and 12 rods are mounted. 
The tubes are heated by nichrome electric heaters. The central tube is inserted through packing, making it possible 
to turn it at any angle without disrupting the operation of the apparatus. 


Twelve nichrome-constant thermocouples, 0.08 mm diameter, are attached to the surface of the central tube. 
The thermocouples are placed in a stainless steel capillary with a diameter of 0.5 mm, The thermocouple junctions 
are calked onto the surface of the tube. The capillaries have been fitted into longitudinal recesses and covered by a 
layer of metal, deposited by means of an EM-Z-A metallizer, after which the surface of the tube was carefully 
cleaned. (The method of making and inserting the thermocouples is described in detail in [11]). 


The temperature of the heat carrier at the inlet and outlet of the heat exchanger is measured by thermocouples. 
Moreover, the temperature of the heat carrier in the individual pile cells at the outlet from the packet is measured 
by thermocouples in capillary covers. 


The tests with mercury and with distilled water were conducted one after the other using the same heat exchan- 
ger. The apparatus for studying heat transfer to water was also made of stainless steel. 


Experimental Procedure and Processing of the Data 

In our tests we measured the temperature field of the surface of the central tube with the heat carrier flowing 
at different velocities. The temperature of the medium at the inlet and the heat flow were held constant. The tube 
was turned from ¢ = 0 to g = 360°, 


For various distances from the entrance,the dependence of the temperature drop from wall to liquid on the re- 
duced angle was determined, For this purpose the data obtained for different central cells were averaged, and the 
angle ¢, equal to 0-360", was transformed to the reduced angles 9, equal to 0-30°. In calculating the temperature 
drops,the mean temperature of the heat carrier in the central cells and in a given cross section was assumed to be the 
temperature of the liquid. A correction for the depth of insertion of the thermocouples was iftroduced on the assump- 
tion that the heat flow distribution over the perimeter of the tube was cosinusoidal. The depth of insertion of the 
thermocouples was ~ 0.2 mm, The average temperature drops over the perimeter were found by graphical integra- 
tion of the relation 9 = f(~), where 9 is the local temperature drop of the wall to liquid (°C). 


In the tests with mercury and water,we used the same central tube I, Control tests with water were conducted 
with a new central tube II. 


The heat-transfer coefficients, averaged over the tube perimeter, were calculated as 


pop aah 
Ql = —= ’ 
iar (11) 
There exists another method for finding the mean heat-transfer coefficient: 
ee cs \ LGN 
IT F pi ee (12) 


P 


But the mean value of the temperature of the heat-exchange surface can be found only in terms of ay: 


1004 


| : a, 
bee == FF \ tak = ti -- de é 
; F ay 
Also investigated in the case of mercury was the heat exchange with a fixed position of the central tube. On 
the basis of these data we also calculated the mean heat-transfer coefficients on the assumption that the temperature 


distribution over the tube perimeter was cosinusoidal. 


In special tests we checked the influence of the unheated case on the temperature conditions of the central 
tube, It was found that heating of the case does not alter the heat exchange in the central cells. This serves as a 
basis for extending the data obtained to piles with a large number of rods, 


The heat flow distribution over the tube perimeter was calculated from Eqs. (6) and (7). 


In processing the experimental data, we took the mean temperature of the heat carrier in the central cells as 
the reference temperature in the similarity criteria. The values of the water and mercury parameters were borrowed 
from {12, 13]. 


Results and Generalization of the Experimental Data 

The limits of variation of the fundamental quantities are given in Table 1. A typical distribution of tempera- 
ture drops is shown in Fig. 2, Such distributions were obtained with mercury flow velocities (w) equal to 1.14, 0.86, 
0.67, and 0.22 m/sec and with water flow velocities of 5.91, 5,80, 5.77, 4.10, 2,22, 2.17, and 1.54 m/sec, The various 
points in Fig. 2 stand for the reading of thermocouples located diametrically opposite one another. The slight diver- 
gence in readings of these thermocouples can be attributed to a small irregularity in heat flowand slight distortion of 
the cell geometry due to allowances in the tube dimensions, It is also evident from Fig. 2 that the temperature of the 
tube surface changes considerably over the perimeter. A similar picture was noted in the water experiments. 


The temperature distribution law for mercury flow- 
ing in the pile is very nearly cosinusoidal. In the case of 
Limites of Variation of the Fundamental Quantities water,more clearly pronounced temperature maxima are 
seen along the lines tangent to the tubes. The tempera- 


AE NBME IE! al 


Value 
Quantity ture distributions, averaged over the cells, are shown in 
for mercury | for water dimensionless form in Fig. 3. With increasing velocity, 
w, m/sec | 0,221,144 1,545.94 the irregularity of the temperature distribution over the 
Lyre 24 SBE Hila perimeter should decrease due to the deeper penetration 
qd, kcal/m?- hr Dil PAO NOE 27,2-103 of turbulence into the narrow parts of the cells. How- 
Asis 0 ee 025 3 saa 9 ever, in the investigated velocity range this effect was 
Re 3800—18600 464016450 only slight. The influence of length was also insignificant. 
Pe 80—470 = 
The value of Tw (A2/A3) averaged over velocity 
and length can be described by the following approximate 
relations: for mercury 
p Ae ; 
Py ~ 0,18 cos 69, (13) 
3 
for water 
T,, 22 ~ 0,124 cos 6p + 0,019 cos 129 + 0,040 cos 18¢ + 0,016 cos 249. (14) 
3 


Figure 4 illustrates the nature of the heat flow distribution over the perimeter of the central tube. The close 
packing is characterized by the fact that radial heat flow along the lines tangent to the rods is completely absent. 
The values of the heat flow at @ = 0°, calculated on the basis of our experiments, are very nearly zero. This is defi- 
nite indication that the results obtained are valid. The heat flow distribution over the perimeter of the central tube 


can be described approximately by the following relations; for mercury 


qiq ~ 1—cos 6g, (15) 
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for water 
gig ~ 1—0,67 cos 6@ — 0,20 cos l2p — 0,07 cos 18 — 0,04 cos 24q. (16) 


In general, the heat flow distributions obtained for mercury and water and basically no different from one another. 


Figure 5 shows the temperature distribution over the length of the pile. It is evident that the heat carrier tem- 
perature in the central cells is appreciably different from the mean temperature over the pile cross section. 


0 30 50 90° 920" 150 TED” ZIP M24 F270 a > S30 Tor 
b Y, degrees 


Fig. 2. Typical distributions of temperature drops over the tube perimeter: a) Mercury, 
w = 1.14 m/sec; b) water, w = 4.10 m/sec. | 


A ANE 2 


Mean Heat-Transfer Coefficients and Nu Numbers, Obtained for Mercury and Water # 


Mercury Water 
(Pr=0,022—0;025) (Pr=3,9—4,9) 
limits of tube a limits of tube | * 
rotation, de cs keal/ Z i Ne i a 2’ : es 
, deg, cal/m?. hr °C rotation, deg, kcal/m?. hr.°C 
0—360 86 2860 0,714 30—210 4650 4910 16,1 
No tum 128 3280 0,82 0—120 6240 7000 23,15 
» » 158 3540 0,89 30—210 5990 5990 19,8 
» » 182 3800 0,96 O0—360 59520 5500 18,5 
» » 204 3850 0,98 30—2:10 11970 9190 30,4 
yo _° 229 4050 1,03 0—210 416300 411680 38,6, 
ay an ne, 1,14 0—90 16450 14850 49,4 
a) ! 4800 ‘bap O0—360 15100 1150 3 
No tum 398 4730 1729 = nes re os: 
a » 439 D000 1,29 — — — = 
0—360 470 5160 1,32 Zs rat = = 
a é ; , 
See Fig. 2 for the angle of rotation notation, Tube I was used in all cases except those marked by b. 
b Tube II was used, 


The mean values of the heat-transfer coefficients over the perimeter a and Nu numbers for the stabilized in- 
terval are given in Table 2, The dependence of Nu on Pe for mercury is shown in Fig. 6. As a result of the relative- 


ly large pile lengths the mean values of the heat transfer coefficients over the length and perimeter differ from the 
stabilized values by no more than 5%, 
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The results obtained for water are compared in Fig. 7 with data calculated from Eq. (1).* This comparison is 
arbitrary with Re = (4-10) - 10°, since Eq. (1) is recommended only for Re > 10‘, The values of Nu calculated from 
Eq. (1) differ from the experimental values approximately by a factor of two. 


\ 
paisa 


Fig. 3. Dependence of the parameter Tw X2/d3 on the reduced angle 9; a) Mercury; b) water. A)l;= 
= 64mm; M)l,=124 mm; @)l3= 184 mm; A)l, = 244 mm; QO); = 304 mm; O)l, = 364 mm. 


The errors in measuring the temperature drops > did not exceed +(0.2-0.3)°C. The error due to geometry, heat 


losses at the ends, and power measurements amounted to ~ +4%, Then for 9 ~ 4°C,the limiting relative errors in 
measuring the mean heat-transfer coefficient and maximum temperature drop from wall to liquid are 


max__ 
6a 442% and ote wey 
Qa 


*In our experiments the correction factor (Prf/ Prw)’’™ barely differed from unity. 
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In processing the results,most of the random errors were eliminated, It may therefore be expected that the systematic 


error 5 (ty *-tf)/9 is less than + 23%, 

Conclusions and Recommendations 
The following procedure is recommended for an approximate thermal analysis of heat-liberating rods and tight 

packing in the pile. The procedure is applicable for Re from 3000 to 20,000. 
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Fig. 5. Temperature distribution over length: a) Mercury, 
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Fig. 6. Experimental dependence Nu= f(Pe) for mer- Fig. 7. Experimental data, obtained for water: O) 
cury: O) Tests without rotation of the central tube; @) tube I; @) tube II; --—--) curve constructed accord - 
tests with rotation of the central tube. ing to Eq, (1). 


The rod parameters should satisfy the approximate simulation conditions (8)-(10): 


w~10 for Pr=0,02—0,03; 
CoC, 130 for Pr=4—5 


(primes indicate natural parameters). 


The heat flow distribution over the rod perimeter is described approximately by the relations (15) for values of 
Pr from 0.02 to 0.03 and (16) for values of Pr from 4 to 5. 


The mean temperature of the rod surface over the perimeter is calculated with the help of the relations 


Nu=f(Pe) and Nu a 
POT (Pry[prors — 1 Re): 
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which are given in Figs. 6 and 7, 


Knowing the heat flow distribution q/q and mean temperature tw, it is easy to determine the local temperature 
of the heat-exchange surface by applying the solution of the heat conductivity problem (6)-(8). 


In tightly packed piles,the dead zones in the vicinity of the rod tangent can cause an irregular deposit of oxides 
and other impurities of the heat carrier on the heat-exchange surface, which leads not only to an increase in the mean 
temperature of the surface,but also to an increase in temperature nonuniformity over the perimeter of the rod. The 
results can therefore be extended only to those cases in which allowance is made for the possibility of contamination 
of the heat exchange surface. 


The authors express their gratitude to A. I. Leipunskii for his lively interest in the paper and for valuable com- 
ments. 
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This paper presents a method of calculating the effectiveness of a partially inserted absorbing, cy- 
lindrical control rod in a reactor without a reflector. The method makes it possible to give a com- 
paratively simple estimate of the distortion of the neutron flux density near the control rod, The 
calculations are made using the one group equation Ag + B’y = 0, with special boundary conditions 
on the surface of the rod. As an example of the use of the method, the values of B” and the corres- 
ponding flux distributions y have been calculated for several positions of the rod in the reactor. 


The effect of a control rod is determined by the amount of reactivity which it counteracts, and depends not 
only on the control rod properties, but on its position in the reactor, 


To get the basic ideas, it is sufficient to look at the simple diagram, shown in Fig. 1. 


In a cylindrical homogeneous reactor of height H and radius R, there 
is a cylindrical cavity of radius a in which a control rod of the same ra- 
dius moves. The control rod is inserted to a depth h(-H/2 Sh =H/ 2). 
The rest of the region, which is not occupied by the control rod, may be 
empty, or filled with the same material as occurs in the homogeneous 
reactor, 


The simplest physical picture of such a reactor configuration is ob- 
tained from the one group approximation. The neutron flux density, 


? = Y(t, z), and the value of the geometric buckling B’, are given by the 
equation 


Fig. 1. Longitudinal section of reactor, 


Ag--B?9 = 0, (1) 


while the conditions at the external extrapolated boundary are 


vei 
p(k, Z=o(1, + )=0, (2) 
and on the surface of the cavity (r = a), the condition holds that 
oP _¢ 
ne K@, (3) 


The values of Rand H/2 are the extrapolated values of these quantities. The derivative on the left side of equation 
(3) is taken along the external normal. The quantity K is a linear operator which depends upon what is inside the 


cavity (r= a). This cavity may be empty, filled with an absorber such as a control rod, or filled with the same ma- 
terial as the rest of the reactor, 


It is difficult to give an exact picture of the operator K, since in any practical case some simple approximate 
value has to be resorted to. Equation (3) holds at the boundary between two media. For an infinite cylindrical cavity, 
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or a completely absorbing control rod in an infinite medium, Eq. (3) can be expressed by the extrapolation length )d, 
which depends on the radius of the control rod, and has a different value, depending upon whether we are dealing 

with an absorbing control rod or a cavity. There are different approximate formulas [1] for calculating AX. One way 
is to consider that the cavity is surrounded by a completely absorbing wall. Then Eq. (3) can be written in the form 


eee Kee —+9. (4) 


If the cavity were filled with the same material as the reactor, then we would write (for R > a), 


2,405 
Ag=> (75 (4a) 


When part of the cavity is filled by an absorbing control rod (inserted to depth h) and the other part of the cavity is 
empty, the approximate expression of Eq. (4) may be used for K. Then 


H 
—6 for —S<2<gh; 


Ko= 
— YQ for Wee Bm 


? 


re) (5) 
Sy 


where 6 and y are the reciprocals of the extrapolation lengths for the cavity and the absorbing control rod, respec- 
tively. If the empty cavity were filled with the same material as the reactor, then instead of the constant 6, we 
would use the constant from Eq. (4a). The function Ky must be continuous at the boundary between regions in which 
we can,with sufficient accuracy,use the asymptotic relation (4). However, it is to be expected that the error will not 
be too great if the asymptotic expression (4) is also used in the transition regions in the vicinity of the point z = h. 

It should be emphasized that discontinuous boundary conditions of this sort are used, for example, in the theory of 
heat conduction. 


From the solution of Eq. (1), with the boundary conditions (2) and (3), it is necessary to determine the minimum 
eigenvalue of B? corresponding with the neutron flux density distribution 9(r, z), Expand the function ¢(r, z) in the 
interval -H/2 =z = H/2 in Fourier's series 


[oe) 


P(r, 2) = > Pr (7) Sn (2) (6) 


k= 
of the complete orthogonal system of functions 
cosa,z for k odd; 
8x (2) = | sin a,zZ for Ak even (a, = r) : (7) 
The cylindrical functions are determined by the following equation 
7.7, B= J, (VB? 2atr) Y,(V B?— of R)— Y, (8? ir) J, (| Bea; R), 
setting 


Zp (r, B) 
Pr (r)= = ®p Z, (a, B?) ’ 


where 9k = 9k (a) represents the as yet unknown constant coefficients for values of k = 1, 2,.... Then for a given 
value of B’ and given values of the coefficients ?,, the function (6) will satisfy conditions (2), 


It remains to satisfy condition (3), which is an operational equation in the space of the function y(a, z) deter- 
mined on the surface of the region. If we represent the operator K on the basis of Eq. (7), then condition (3) will be 
equivalent to the system of linear equations 


H Zh(a, B?) ae a 
ae k= >} AQ (k= Les sa) (8) 


tt 


IO HEL 


g i K as determined by (7). 
where A Ge Loe — Lp aay he | eee and Kyy are matrix elements of the operator y (7) 


Thus the problem is reduced to the solution of the infinite system of equations (8), which determines the eigenvalue 
B? and the corresponding eigenvector 9k. 


We now use the approximate expression for the operator K as given in equation (5). For the elements of K,we 
obtain the following expression 


H 
h a 

— Kj, a 6 \ dzgy, (2) 8) (2) ae ni \ dzg,, (z) £1 (2). (9) 
H v 


Mee 


The minimum eigenvalue of B? according to Eqs. (8) and (9) is determined from the equations 


iE 
Zi(a, BY) =y0,(a, BY) (k= —> (10) 


a 


and 


Zi (a, B?)=6Z, (a, B?) (h=+). (11) 


Equation (10) holds incase of a completely inserted absorbing control rod, and Eq. (11) holds for an empty region. 
The flux density distribution in both cases is given exactly to an arbitrary constant multiplier, by the function 
Z, (rt, B’) cos a4Z. 


The determination of the eigenvalue and the corresponding eigenvector from the infinite system of equations 
(8) is almost as difficult as the solution of the original problem. In practical cases, however, the following method 
of solution is possible. Set yk = 0 for k > N, and replace Eq. (8) by the finite system of homogeneous equations 


H N 

al 

= Fe a= Dy Kids (lest Faecal (12) 
i 


This finite system determines an approximate value of B’ and the corresponding values of the Fourier coefficients 9, 
(where 1=k =N), from which Eq. (6) gives the approximate distribution of the neutron flux density. 


The problem is considerably simplified by the fact that the value of B’ lies somewhere within the interval 
Bins B’< oor the end points of which are determined by Eqs. (10) and (11). 


a lnieael 


Geometric Buckling B” and Eigenvector ¢xK for a Control Rod Inserted to Half Depth (h = 0) 


N |B2.103lqy ¢2 | {3 | G4 | 5 6 | Q7 | Ps Qo | P10 | Pil 
| 
4 0, 2078] 1 |—0,865 |—0,308 | —0,0933 = == eee = aoe a =e 
9D |0,2064| 1 |—0,844 |—0, 216} —0,231 | —0,184 = 22 = ota == Se 
6 }0, 2064) 1 |—0,838 |—0, 209 | —0,224 | —0,160 | —0,0374 -- — — — oe 
7 10,2054! 4 |—-0,824 |—0,173 | —0,242 | —0,121 | —0,120 —0,148 — — — — 
8 |0, 2054] 1 |—0,822 |—0,171| —0,240 | —0,118 —0,117 | —0,105 | —=0,0209 — — — 
9 0, 2049 4 |—0,814 |—0,152} —0,245 | —0,0996] —0,132 —0,0835) —0,0767} —0,0851 —_— — 
i ees —0,812 |—0,151 | —0,245 | —0,0983] —0,134 —0,0813] —0,0747 —0,0774 —0,0139 — 
yo 


—0,806 |—0, 139 | —0,247 | —0,0876] —0,137 | —0,0700] —0,0866 ~—0,, 0628) —0,0547} 0,0656 


The solution of a finite number of equations has been carried out for a reactor of radius R = 215 cm and height 


415 cm with a cavity of radius a = 2.5 cm and the reciprocal extrapolation lengths 6 = 1.726 - 1074 cm™}, y = 1.1629 


7 
cm’ for three values of h; namely, h = H/4, h = 0, and h = -H/4, i.e., with the control rod inserted one quarter, 
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one half, and three quarters of the length of the active zone. For a reactor without reflector, having these same di- 
mensions, the buckling is 


B= (# by + (=e) = 0s18 2454 028, cmt 


Since the radius of the cavity, a, is small, and the value of 5 is very near zero, the value of B? from Eq. (10) is prac- 
tically equal to the value of B’ for a bare reactor without a cavity; i.e., Bonin ~ 0.1824: 10°? cm72, The value of 
Beer from Eq. (11), corresponding to a completely inserted control rod is 0.2244 - 107? cm ”, 


Now it can be shown that of the functions Z’, (a, B’)/ Zk (a, B’), only the first (k = 1) depends upon B? in an 
essential fashion, and the rest are almost independent. Consequently, if in formula (12) we leave out the first equa- 
tion, and solve the system of the remaining N-1 equations by any numerical method for any value of B’ in the in- 
terval Bifbo abc Bea) and 9, = 1, we will get an approximate value of the eigenvector ¢ of the finite sys- 
tem (12), After substituting the gy, into the equations of the system, we can calculate the corresponding eigenvalue 
of B’, starting from which the values of ~, and B’ can be improved by further iteraction of the solution. For calcu- 
lations of this sort it is convenient to have a table prepared ahead of time for the ratio Z'; (a, B’)/ Z,(a, B’) as a 
function of B? in the interval (Bonin Brash 


Table 1 gives the values of B’ and the corresponding eigenvectors obtained by the solution of a finite number 
of equations for the case h = 0. It is clear from the values given that with this method of calculation, B converges 
fairly rapidly with increasing N. The flux density distribution g(a, z) on the surface of the cavity is given for four 
values of N in Fig. 2 a andb. The calculational results for control rods inserted one quarter and three quarters of the 
depth (h = + H/4) are given in Tables 2 and 3. The curves of the flux density distribution g(a, z) on the surface of 


WAV be 2 


Geometric Buckling B’ and Eigenvector Yy for a Control Rod Inserted to One Fourth Depth (h = H/4) 


| N |B2.103| gy | Qe | 3 Pa | 5 | 6 | Q7 Ps | Po | 910 
4 |0,1889 1 =), Piks} 0,180 0,0895 | — — — — = a 
6 |0,1881 1 —(),249 (),184 0, 0869 0,0127 | 0,0775 = = 22 Ae 
8 10,1874 | (0), 1s 0,218 0,0793 |—0,0108 | 0,0274 —0,114 0,0126 — — 
10. 0, 1868 { == (Oe Sltei7/ 0,218 0,0842 |—0,0200 | 0,0221 —0,121 |—0,04117 |—0,0176 | —0,0378 
TABLE 3 


Gometric Buckling B’ and Eigenvector yx for a Control Rod Inserted to Three Fourths Depth (h = ~H/4) 


N |B2.103| gy | G2 (3 pa | 5 | 6 | —7 | 8 | 9 | 16 
Aue 24 1 Sy AO) || SOR AG 0,380 = = — = = _ 
6 10,2224 1 == ()) 600m e=—O nee 0,642 0,448 —(), 220 — = = = 
8 10,2222 4 ==(0 6061) 90, (08 0,677 0,483 —0,251 |—0,0983 | —0,0538 = —- 
OM O22 241 1 —(),687 | ==0), 769 0,670 0,443 —0,179 |—0,00443/ —0,155 | —0,165 0,105 


the cavity for these conditions are given in Figures 3 and 4. It is clear from the curves that the distortion of the 
neutron flux along the z axis can be quite large as compared with the original cosine distribution on the surface of 


the control rod. 
Since the function Z,(r, B’)/ Z(a, B’) rapidly approaches zero with increasing r for k > 1, and the more ra- 


pidly the greater the value of k, the higher harmonic terms in series (6) diminish relatively quickly. At a certain 
distance from the cavity, the distribution along the height is practically determined by the first harmonic, the cosine. 


This condition is illustrated in Fig. 5 (h = 0, N = 4). 
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Fi 6 shows the geometric buckling B’ as a function of the depth of insertion of the control rod, hr. 
igure 


Pp y i i i y i ri al con 
| 1 * . ° ) 


Fig. 2, Flux density distribution on the surface of the cavity (h = 0); a) 
ist Opies eae Te, 4s b) ————._ N= 10, —- — — — Nive 


The one group approximation is used in this paper. A generalization of the method to the two group case would 
not present any fundamental difficulties, but one would have to contend with the difficulties which are peculiar to 


the two group approximation, such as determining the proper boundary conditions for the fact neutron group at the 
control rod surface. 
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Ad 
2 a 3 


Fig. 3. Flux density distribution on the surface of the cavity (h = H/4): 
)N=10; ————— )N=6. 


)N= 10; ————— ViNi=IGs 


ail 
4 


Fig. 5. Flux density distribution at different distances from the cavity, with the 
control rod inserted to half depth (h = 0). Distance from the cavity (in centime- 
ters); 1) 2.5; 2) 10; 3) 20; 4) 30; 5) 100; 6) 150. 
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Fig. 6. Geometric buckling as a function of 
depth of insertion of control rod. 


The results of the calculation of controlrod effectiveness, which we have arrived at in the present paper, are to 
be compared with the calculations given in reference [2]. Although these results do not lend themselves to direct com- 
parison, the curve shown in Fig. 6 is in complete qualitative agreement with the curves of reference [2]. 


We wish to thank R. Zezula, I. Svatosh, and M. Prazhskaya for carrying out the numerical calculations. Parti- 
cularly great initiative was shown by R. Zezula, who also calculated the tables of the function Zz (a, B’y/ Zh (a, B’) 
on the electronic BESM setup in the Calculational Center of the Academy of Sciences of the USSR. 
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INVESTIGATION OF THE ISOTOPIC COMPOSITION OF URANIUM 
IN RARE-EARTH MINERALS 
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December, 1960 
Original article submitted February 24, 1960 


The isotopic composition of the uranium of rare-earth minerals (orthite, xenotime,and gado- 
linite) was investigated for the detection of possible traces of Cm™’ in nature,by means of an 
ionization chamber with a multichannel analyzer of the pulse amplitudes. The results of the 
measurements are discussed. A method is given for the accurate, calculation of distortions of 
the spectrum of @ particles in ionization chambers with electron collimation. It is shown 
that during the investigation of samples of low intensity (0.2-1 mg of uranium),the relative 
content of U*> /U™® isotopes can be measured by an a-spectrometric method with an error 
up to ~ 1%, 


Introduction 

An investigation of Cm™" obtained in a reactor by irradiation of plutonium [1] made it possible to determine 
the lower limit of the half life of « Jecay, equalto4-10' years. The absence in Cm" of types of decay, other than 
spontaneous fission (Tf > 6 - 10° years), indicates that it may exist in nature. The discovery of the highly anomalous 


content of U® in magnetite by the authors of [2] is also in favor of this. 


In its geochemical properties, curium is an analog of the rare-earth elements and, if it does exist in nature, 
would most probably be found in ancient rare-earth minerals. 


A rough assessment of the amount of Cm" which can exist in nature at the moment of formation of minerals 
indicates that its direct determination in rare-earth minerals is only possible if its half life is not less than ~2- 10° 
years. But if Cm” existed in nature and decayed during geological time, then, being converted to wee by the 
following degradation 


a B- a Biz a 
47 2 ‘ a 
Cm?4? —> Pus ate Am arnt > Np??? Saas Put? meee aiercre Se 
h 600 yr 2, ay 24400 yr fishe tO? yt 


it must have changed the isotopic composition of uranium (u”> /u™® ratio), the change being particularly apprecia- 
ble in the case of rare-earth minerals having a slight uranium content. 


In the present work an investigation was made of the isotopic composition of uranium in ancient rare-earth 
minerals, xenotime, orthite,and gadolinite, the age of which is ~ 2 billion years, i.e., about the same as the geo- 
logic age of the earth, Samples of these minerals were obtained from the Mineralogical Museum, Academy of Sci- 
ences of the USSR. The uranium content in xenotime, orthite, and gadolinite was 0.8, 0.1, and 0,06%,respectively. 


Radiochemical Separation of Uranium from Minerals 

For the a@-spectrometric investigation, uranium was separated from the minerals by the following method. Af- 
ter the acid decomposition of the investigated minerals (HC1 for orthite and gadolinite, and HSO, for xenotime),the 
solutions obtained were evaporated to dryness, nitric acid being added. Nitric acid and magnesium nitrate were then 
added to the residue until the concentration was 0.5 N and 2,5 M,respectively (a small amount of iron nitrate was 
also added to the solution obtained from xenotime, to combine the phosphate and sulfate ions), The uranium was 
then extracted with three equal amounts of ether. The ethereal layer was separated and evaporated to dryness in the 
presence of dilute hydrochloric acid. The residue was dissolved in 8N hydrochloric acid. After traces of iron had 
been removed by extraction with amyl acetate, the hydrochloric acid solution was passed through a chromatographic 
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column filled with Dowex-1 anionite; the latter was washed with 8N hydrochloric acid, while the uranium was eluted 
with 0.5 N hydrochloric acid. The solution obtained was evaporated to dryness and chromatographic purification was 
repeated. The purified uranium was precipitated electrolytically on stainless steel discs (diameter 70 mm) by a meth- 
od described in detail in [3]. The hydrochloric acid solution of uranium was evaporated almost to dryness, and the 
residue was dissolved in a 1% solution of ammonium oxalate. This solution was then acidified with 6N nitric acid. 
Electrolysis was carried out at a temperature of 80°C and a cathode current density of 100-150 ma/ cm”, During elec- 
trolysis, several drops of 6N nitric acid were added periodically to the solution. After 15 min the solution was made 
alkaline with several drops of 6N ammonia, and electrolysis was continued an additional10 min, The solution was 
mixed during the whole period of electrolysis. Under these conditions, uranium was 80-90% precipitated. 


The preparations obtained by this method were perfectly suitable for &-spectrometric measurements. 


Measurement Method 

The relative content of U2 and U”* isotopes in the minerals was determined by the activity of these iso- 
topes. Since the amount of uranium in the samples investigated was small (0.25-1 mg), an apparatus with maximum 
possible illuminating power — a pulse-counting ion chamber with a grid — was used. Figure 1(dashed curve) shows the 
a spectrum of the uranium isotopes, obtained in such a chamber. But the u”* content cannot be determined exactly 
by means of this method because of the considerable "tail" from 
the U4, This tail is due to the absorption in the layer of the source 
of part of the energy of the a particles escaping parallel to the elec- 
trodes. The quality of the spectrum can be improved by using a 
mechanical collimator. In the first place, however, existing col- 
limators do not allow sufficient improvement of the quality of the 
spectrum and, secondly, markedly reduce the counting rate. The 
use of electron collimation gives much better results. Various meth- 
ods of introducing electron collimation are described in [4-6]. They 
are all based on the use of the relation between the pulse height at 
the high-voltage electrode and the angle of emergence of the a 
particle (height method of collimation) or the time lag of the pulse 
at the collecting electrode with respect to the pulse at the high-vol- 
tage electrode (time method of collimation), Use of both methods 
gives a substantial improvement in the quality of the spectrum with 
a relatively small decrease in the counting rate. In Fig. 1 the spec- 
trum obtained by means of time collimation is represented by the 
solid curve. The counting rate here was 100% less. From a com- 
parison of the spectra,it is evident that as a result of the use of col- 


Number of counts in the channel 


0 0 F v5 se = ator the qagiiry of the spectrum is considerably improved, but 
Numberor the channel that during this process the relative intensity of the lines is distorted, 
Fig. 1. Alpha spectrum of a natural mixture whichis a shortcoming of electron collimation, It is found, how- 
of uranium isotopes, obtained during measure - ever, that this distortion can be taken into account with high accu- 
ments with (solid curve) and without (dotted tracy. To calculate the corresponding correction we will examine 
curve) collimation. in greater detail the time method of collimation, used in the pre- 
sent work. 


Figure 2shows the diagrammatic layout of the ionization chamber. The o-particle source is located on the high- 
voltage electrode. The emitted particle ionizes the gas with which the chamber is filled (A + 0.5% CH,). Under 
the effect of the applied electric field the electrons begin to move toward the grid, which leads to the immediate 
appearance of a pulse at the high-voltage electrode, irrespective of the direction of escape of the a particles. A 
pulse appears at the collecting electrode only after the electrons nearest the grid have reached the latter. The time 
lag tlag is evidently independent of the orientation of the track and can be expressed by the formula 


d— Roos 
tag= ’ (1) 


where R is the o-particle path, w is the rate of drift of the electrons; d is the distance between the high-voltage elec- 


trode and the grid; y is the angle between the direction of emergence of the a particle and the perpendicular to the 
surface of the electrodes. 
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It can be shown that the distribution of the number of pulses with respect to the time lag must be uniform (with 
isotopic distribution of the a particles), as is shown in Fig. 3 for two groups of o particles of different energies. The 
a particles emerging perpendicularly to the surface of the electrodes correspond to pulses with minimum lag. 


The method of time collimation is based on the fact that only those pulses for which tlag is less than a certain 
value of t’ are recorded. The pulses from a particles emerging at low angles are not recorded. The degree of col- 
limation can be varied at will by varying the value t’. 


Let us determine the variation of the ratio of the line strength 
as a result of the introduction of collimation. The degree of col- 
y  limation can be conveniently characterized by the value f, deter- 
mined from the ratio 


= 
x f ' 
wo wy —— tmax— t . (2) 
ass Sess Eee See eo Hence, it is evident that the proportion p of recorded pulses of a 
5 particles of a given line is determined by the relationship 
oO 
Fi 2 Di . . . . . & No a ; i R : 
a * Sia SSE: son aon where Ng is the line strength before the introduction of collimation; 
oe a) Den voltage electtoceis2) SUGk = Nis the strength of the same line after the introduction of collima- 


3) collecting electrode. rion 


From equation (3) it is seen that the variation of the line 
strength depends on the path (and, therefore, on the energy) of the 
a particles. 


= 


Let there be two lines with strengths NY and N}. In the experi- 


2 ment with collimation, the ratio Q = N,/Ng is measured, but it is 
a. 0 _ p10 74,0 
necessary to measure the value Q: = Ny/Np. 
Se 
fo) 
= Employing equation (3) we obtain 
te) 
ag 
=) 
Z rpc ie (a 
OSS Ni a Ny Ry _— Ry 4 
SORE ND GN ip Toa raeai aatmted g a 
ie iis 
Fig. 3, Distribution (theoretical) of pulses where A is the correction factor. 
i i f f 
oS Dod epee us pe yoeges te ees g If the degree of collimation (i.e., the value f) is known, the 
a particles of different energies and intens - } ; Gees : 
: true ratio of the line strengths Q with respect to the observed Q can 
Bes: be found by means of equation (4). In consequence, the whole prob- 


lem amounts to the accuracy with which the degree of collimation 
can be determined. 


During the investigation of the isotopic composition of uranium ores,the following method of determining the 
correction factor was found to be simple and reliable. In the given case there are three lines in the spectrum, two of 
these (U* and U*) being fairly strong, and therefore suitable for measurements. Subsequently, all the values relat- 
ing to U™, u™> and U™* were given the indices 234, 235,and 238, respectively, and, in addition, the following sym- 


bols were used for the ratio of the line strengths: 


Qe _ N2ss Q __ Nose Qe _ Nase a) ee Noga 
234°> = ’ 234) ’ = ’ 9s. — on 
N34 4 Noga a IN Bs 5 Nogs 


For the investigated ratio of the line strengths of u? and U2, equation (4) is converted as follows: 


Qoss = Qess4eas = Ve35 a8 Qoss (1 + Hag Bese = Qogs (1 + Poss), (5) 


Qos Qo38 
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where 


— Qons—Qens at) F 
Poss = ied pee Rezp Ros (oa) 


Similarly, for the u? and U** lines we have 


Qos = Qosa4eaa = Vesa (1+ Posa), (6) 
where 
vat f 4 4 
Poss 7 4 t ( Rog Rosa ) : (6a) 
Roge 


The values of QS, and Qoy are determined by experiment, which makes it possible to determine Aggq and, therefore, 
Pog. Further, by employing equations (5a) and (6a),it is easy to establish the relation between the Pogs and Pog4: 


4 4 
_ Rog Rass __ Rogs—Rose Rosa p  — bP ne (7) 
ae 4 eee! Pasa = Roga— Roaa A235 a ; 
Roz Rega 
where 
b a —tes5— Ress Rose (8) 


’ 
Rosa — Rags Ross 


Since the paths are only present in equation (8) in the form of ratios, the coefficient is determined only by the energy 
of the a particles, and is independent of such values as pressure, temperature, etc. During the determination of the 
numerical value of the coefficient,the differences of the paths can be replaced by the difference of the energies. The 
ratio Ry34/ R35, determined by the path — energy curve [7], is equal to 1.135. The following values are adopted for 
the energies of the a particles: Eggg = 4.195 Mev, Egg, = 4.767 Mev, Eo35 = 4.390 Mev. Substituting these values in 
equation (8), we obtain b = 0.39. 


TABLE 1 


Results of Measurements without Collimation 


Amount of |Measurement 


0 0 0 
Sample uranium, mg |time , hours Nesa | Nass | es 
| 
Standard, urnayl nitrate ........ | 1,50 4,0 118 200 118 750 1,005 
No, 1, Omthite Weak ae eek 0,65 4,5 42, 280 41 800 0,989 
PCR COOLING cue austere Sess aa 5's 0,42 8,0 47 030 47 510 1,040 
4,0 40 150 37 800 0,944 


Re: Or tauOliilic se Ae eee 0,18 4. 


Instead of the ratios of the energy differences, the ratio of the distances between the lines in the spectrum in 
Fig. 1 may be used. The same value is obtained for the coefficient b in this case. 


By means of equations (6) and (7) we determine the correction factor for Ao3s: 


Aas = (1 + bP a4) = [4 + b (Agay — 1)] = [4 + 0,39 (Aas — 4)]. (9) 


Thus,for the investigated ratio Q2ss we obtain the final equation 


Qos = Qos5 [1 + 9,39 (A,3, — dy): (10) 
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It can easily be shown that the relative error 6 Ags is related to 5b by the following equation: 


_ 9b [b (Aog4—1)] + 6493440950 
OA gan =e ee Ong h2340 
as i+b (dg) ae 


For the degree of collimation, Aj * 1.20, actually employed, 
6Ag35 y= 0,165 -- 0,06 Aga4. a 


The coefficient b is determined with an accuracy up to 5%, This means that the part of the error 5 Aga, due 
to the inaccuracy of the determination of b, does not exceed 0.5%, It may be noted that this error is systematic and 
occurs only during absolute measurements. In the case where the source investigated is compared with the standard 


and they are both investigated under approximately equal collimation, this error is still further reduced and may be 
neglected. Hence, 


8 Asgs © 050Aoa,. 


The value of 6A234, determined only by statistics, can amount to a fairly small value (0.5-1%), depending on the ac- 
curacy required. 


From the above it follows that to determine the cor- 
rection factor Aggs,it is sufficient to know the ratio of the 
line intensities of U?* and U? before and after the intro- 
duction of collimation; i.e., the values Opes and Qog,. 


To determine the value of the ratio Qoy,it is neces- 
sary to carry out a further experiment; it is obtained direct- 
ly from the experiment with collimation for the determi- 
nation of the U® content. Only the experiment without 
collimation for the determination of the value of the ratio 
Ghen is an additional experiment, the time spent on this 
being roughly five times more than for the principal meas- 
urement. The fact that the correction factor Ag3s is deter- 
mined directly from the principal experiment is all the 
more useful because,in this way, possible variations in the 
degree of collimation during the measurement period are 
unimportant. 


1500 


1200 


900 


600 


Number of counts in the channel 


Thus, the alteration in the ratio of the line intensi- 
ties, caused by the introduction of electron collimation, 
may easily be allowed for by an error not exceeding 0,2- 
0.5%, 


300 


Results of the Measurements 

Number of the channel To check the method and the agreement of the 
Fig. 4. Alpha spectrum of the investigated samples of measurement results, the y/ U8 ratio was determined 
uranium: ) standard sample;(X)sample No. 1; O) for uranium samples obtained from commercial uranyl 
sample No. 2; A) sample No. 3. nitrate and uranium minerals, The measurements showed 

close agreement of the results, The mean ratio ~ 1/138 

was taken as the unit. Subsequently, all measurements were carried out with respect to one sample, taken as the 
standard. The value q = Oss standard/ Qdss sample Was determined, 


The value of q was determined as follows, The ratio Qbe4 was first obtained. To increase accuracy, this meas- 
urement was carried out with a mechanical collimator with a low degree of collimation (aluminum disc, 2 mm thick, 
with an orifice of diameter 4mm). The main measurement was then carried out with electron collimation. The 
degree of collimation was established in such a way that the counting rate was reduced by roughly 100%, The ratios 
Qos, and Qoss were determined from this experiment; then Qo35 was calculated by means of equation (10). 
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The measurement obtained for the background showed that it could be neglected. The low value of the back- 
ground was obtained because the pulses from the high-voltage and collecting electrodes, coincident in time, were 
recorded, The results of the measurements and calculations for the standard and the three investigated samples are 
given in Tables 1 and 2. Figure 4 shows the® spectrum of the investigated uranium samples. The sPeene Rae a com- 
mon height in the maximum of the u* line. As may be seen from Fig. 4, although the tail of the U ie we mar- 
kedly reduced as a result of the introduction of collimation,it is nevertheless 2-3% Ngs. The complete similarity of 
the investigated spectra (samples No, 1 and 2 and the standard) makes it possible to affirm that,in the geal 
this proportion is equivalent to an error of not more than 0.5%. Bearing this in mind, we conclude that the U™™ con- 
tent in samples No. 1 and 2 was determined with an error of less than 2%, As regards sample No. 3, the measurement 
error may be expected to be somewhat impaired by the presence of a small amount (traces) of Th”, But in this case 
the separation of the uranium isotopes was,to some extent, even better,which is explained by the more careful pro- 
cessing as a result of the deficiency of the mineral of the laminated sample. In this case, too, the measurement error 
of the U*>/U”*8 ratio is, therefore, ~2%. 


As may be seen from Table 2, a somewhat higher isotopic u?>/u?™ ratio (q 1.046 + 0.02) is observed in 
gadolinite. This indicates that the uranium of the given rare-earth mineral contains an excess of U™> and equals 
~ 0.03% of the total uranium. As already indicated, a possible explanation of the anomaly discovered in the isotopic 
composition of uranium may be the existence of curium in nature, 


TABLE 2 


Results of Measurements with Collimation 


Meas - 
Sample neat Noga Nog8 Ness Qesa A234 Qe35 or q 

time, 
hours 

Standard, uranyl 

Nitrate: cee ceeue ele. 30 514 100 | 408070 18 124 0,795 1,264 22,50 24,87 1 

INGey, th CNRS. 5 G6 6 20 164 660 | 138 080 5 864 0,839 1,179 23,00 25,20 0,987 

No, 2, xenotime... 30 157 880 } 130590 5 706 0,827 124 22,89 24,86 4 

INO, Sh. gadolinite, Size 70 448 310 | 117670 5 310 0,793 1,187 22,16 23,78 4,046 


If it is borne in mind that the age of the gadolinite investigated is ~ 2 billion years and the half life of Cm™? 
is 24-10! years, then, neglecting the time required for conversion of all the curium to U™®, the amount of curium 
present in the investigated mineral at the period of its formation can be roughly assessed. The initial content of 
Cm™" in the given sample of gadolinite with a mean uranium content of 0.06%,is 2 107%, 


In the near future we plan to use the new method to improve the accuracy of the data obtained and carry out 
an investigation of other minerals. 


In conclusion, the authors with to express their thanks to A. P. Komar and V. I. Baranov for the interest they 
have shown in this work. 
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METHODS OF DIRECTIONAL DETECTION OF GAMMA RADIATION 


View bee B Onvaln 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 483-487, 
December, 1960 
Original article submitted February 22, 1960 


An analysis is given of the principal characteristics and parameters of directional radiation de- 
tectors. Particular attention is paid to directional scintillation counters of the compensated type, 
and it is shown that they possess much better parameters than analogous instruments using gas- 
discharge counters. A description is given of various compensating systems which may be useful 
in connection with the detection and measurement of radioactive elements. 


Directional recording of y radiation is being increasingly used at the present time.. This type of detection is 
employed whenever it is necessary to determine the position of the source of radiation or to measure the intensity in 
a given direction in the presence of other interfering radiation [7, 9, 10]. 


The distinctive features of a directional radiation counter are best shown by a polar diagram which gives for 
each direction (y, 9) the ratio of the y-ray sensitivity q(y, 9) of the detector to the y-ray sensitivity in the direc- 
tion y = 0, 9 = 0; ice., 


_ 9(@, §) 
n(Q, DS rr (1) 


where ¢ and 0 are the two angles defining the direction of incidence of the y rays. 


The polar diagram is usually given either in the 9 plane or in the gy plane. The angular resolution of sucha 
system is usually described by the angular width of the polar diagram o (Fig. 1),which is defined as the angle be- 
tween the two directions corresponding to one-half of the maximum sensitivity. 


Directional radiation detectors can be used to measure weak signals against a strong background, owing to the 
improved signal-to-background ratio. In fact, the root-mean-square relative error in the counting rate np in the pre- 
sence of background nb > ny is given by 


5 not 2np/ Re Ant 2n, (2) 


where T is the time of measurement,and Rg = no/ ng is the directivity. 


The directional y-ray sensitivity characteristics can be obtained in various ways. The main types of radiation 
detectors which have been used in this connection are special gas-discharge counters with bimetallic and ribbed cath- 
odes and insulating beads on the anode [1, 2], composite radiation detectors incorporating coincidence and anticoin- 
cidence circuits [2, 3], radiation detectors with screens and collimators [4, 5], and compensated systems [6, 7, 10]. 


Special types of gas-discharge counters and composite radiation detectors have a rather low directivity and 
efficiency, and their characteristics are very dependent on the energy of the radiation. They have,therefore not been 


used in practice. 


The most widely used directional detectors are those of the screened type, because of their simplicity of con- 
struction and high directivity. The thickness of the screen, its form, and the dimensions of the aperture determine 
the angular resolution of the radiation detector. As a rule, the screens are made of lead and their thickness is suffi- 
cient for the virtually complete absorption of radiation incident upon them. Significant disadvantages of screened 
detectors are theirlarge weight and the dependence of the directivity on the energy of the radiation. Thus, for .exam- 
ple, the screen described by Takahashi et al, at the Second International Conference on the Peaceful Uses of Atomic 
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Energy in Geneva [7],had 4 cm thick walls, an angular resolution of 10°,and a weight of about 160 kg. Even with this 
screen thickness, radiations with energies in excess of 1 Mev were only partly absorbed. 


The effect of unwanted radiation can be completely excluded in directional systems of the compensated type, 
in which the fraction of the unwanted radiation passing through the screen is compensated with the aid of an addition- 
al radiation detector. In this method, signals from both detectors are recorded in separate channels (Fig. 2a) and are 
fed into a differential counting device. With the radiation incident in a direction ¢, the latter device produces the 


difference signal 


Ang = Nig — Nag, (3) 
where nig = mqigN, neg = mqggN; N is the number of y rays of given energy reaching the radiation detector per- 
unit time and m is a constant. The radiation detector having a large sensitivity in the direction y = 0, for which Ang 
is a maximum, will be referred to as the main detector. A directional sensitivity characteristic can be obtained with 
the aid of two types of compensating systems. In the first type (Fig. 2b) the main detector is screened, while in the 
second type (Fig. 2c) the compensating detector is screened. In order to compensate the fraction K of the unwanted 
radiation passing through the screen and recorded by the main detector, the sensitivity of the compensating detector 
in the system of the first type (Fig. 2b) should be of the form 


q2= 1K. (4) 
In this case,a parallel beam of y rays in the principal direction (y = 0) gives rise to the difference signal 
Any = Nyy — Ngo = No (1 — K?). (5) 
In systems of the second type (Fig. 2c), the two channels are adjusted to the same sensitivity and hence the dif- 
ference signal recorded for a parallel beam of y rays incident along the principal direction is given by 


Any = Nyy — Nog = Ny (1 —K). (6) 


Comparison of Eq. (5) with Eq. (6) shows that the sensitivity of directional detectors of the first type is higher 
by a factor of 1+ K, The directional properties of a compensated system can be represented graphically as the dif- 
ference between the directional characteristics of the main and the compensating detectors. 


180° 


a) 


90° 


ZI 


Aus 


Y, 
Fig. 1. The directional sensitivity of a compensated system of the first type 
in polar and Cartesian coordinates, 


Figure 1 shows the directional characteristic for the compensated system of the first type (cf.Fig. 2b). 


A disadvantage of both types of compensated systems is the increase in the statistical fluctuations at the differ- 
ential counting unit, which are determined by the total error in both channels: 


o=V)Vo'+o3, (7) 
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where 6, and 0, are the root-mean-square errors in the main and the compensating channels, 


Let us consider the characteristics of the two types of compensated systems in an isotropic y -radiation field. 
The number of pulses An recorded per unit time by the differential counting unit is given by 


AN = \ n(g, 9) dQ, (8) 
Q 
where d& = sin 6 dO dg is an elementary solid angle, n(¢, @) is the number of pulses recorded by the system in the 
direction (¥, 9) per unit time,and Q is the solid angle within which n(¢, 9) differs from zero. 


The function n(¢y, 9) can either be calculated or de- 
termined experimentally. When it is unknown, and only 
Main Main channel I, ; the total counting rate An and the counting rate Ang for 
detector radiation incident in the principal direction only are given, 
it is convenient to use the effective solid angle ®. This 
angle characterizes the directional properties of an equi- 
valent system with a constant sensitivity within this angle 
and is defined by An; i.e., 


An — An,4@,where @D = =a = 
0 
Ree ck (9) 


2 


The polar diagram for this case is in the form of a spheri- 
cal sector with a radius equal to Anp and a volume nume- 
rically equal to An. 


Using the concept of the effective solid angle in the 
case of the compensated system of the first type, the num- 
ber of y rays recorded in the main channel per second is 
given by 


Fig. 2, Functional scheme and the two main types of 
compensated radiation detectors with directional prop- 


erties. ny=nO+nk (1—9), 


where n is the number of y rays recorded per second by the main detector with an isotropic sensitivity characteristic. 


In the compensating channel, the number of counts per second is given by ng = K [Kn® + n(1-®)]. The root- 
mean-square error is then given by 


HV 24+ B=) 2 [oO — Ky 42K), (10) 


where T is the time of measurement in seconds. 


The differential counting unit will record An = ny-nz = n® (1-K4 pulses/sec . It follows that the relative root- 
mean-square error is given by 


A consideration of Eq. (11) shows that the error rapidly increases for screen transmission coefficients greater 
than 0.5 and for low effective solid angles. This is equivalent to a large unwanted radiation intensity, or a reduced 


fraction of radiation recorded by the compensated system (Fig. 3). 


In the case of a narrow beam of y rays incident along the principal direction, it may be assumed that ® = 1, 
In this case Eq. (11) is considerably simplified and the error is given by 


5 1 yee 
1" 4— K? nT —~ 
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In the case of the compensated system of the second type, the fluctuations will be large,since the unwanted ra- 
diation is recorded by both detectors without attenuation and the sensitivity in the principal direction will be lower 


by a factor of 1 + K. 


In fact, using the same notation as in the previous case, we find that 


5 or) 1 


Z An N=TK 


pai tl n.nkKO{n(i—®) _ , /2n—nO(i—K) . 
i yj = spe Vrs T = V ik Oo 


Ve ase (13) 


Even in the most favorable case,when ® = 1, the error due to statistical fluctuations in the compensated system 
of the second type will be greater by a factor of [(1 + K) V(K + 1)/(K? + 1)]. 


In view of this, systems of the second type are recommended when it is necessary to reduce the weight of the 
detector to a minimum and in the case when it is necessary to record local sources against relatively low unwanted 
background. These systems are advantageous even in cases where the detectors have a large background np due to 
cosmic radiation and radioactive contamination. When the sensitivity of the main and the compensating channels 
are the same, this background has no effect on the differential count, while in systems of the first type there is a con- 


stant component equal to nb (1-K). 
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Fig. 3. Dependence of the relative root mean square 
error on the screen transmission coefficient and the 
effective solid angle for a compensated system: ---) 
nT = 100; ——) nT = 1000. 


Gas-discharge counters or scintillation counters can 
be used in compensated directional detectors. Analysis shows 
that in the large majority of cases,scintillation counters have 
important advantages owing to their high efficiency, good 
time resolution, small phosphor dimensions,and the possibil- 
ity of selection of spectral sensitivity characteristics. 


In experimental studies of the efficiency of gas-dis- 
charge and scintillation counters incorporating various kinds 
of screens [8], it was found that in the case of the scintilla- 


tion counters it is possible to select filters and working con- 


ditions so that Eq. (4) is satisfied in a sufficiently wide ener- 
gy range. 


Figure 4 shows a number of compensating systems and 
their directional characteristics. The scheme with two ra- 
diation detectors having equal sensitivities and a screen Ky 
between them has a two-lobe characteristic (Fig. 4a). The 
sensitivity in the reverse direction (y = 180°) can be reduced 
by an additional filter Ky. Without this filter the system can 
be used in the absence of radiation in the directions 90° < 
5” @.<270:, 


The detector shown in Fig. 4b has a complex direc- 
tional characteristic ensuring the same sensitivity in opposite 
directions along the x axis and one of the z directions, As 
a result, a Zonal characteristic is obtained which is conven- 
ient in the location of moving local sources of radiation. The 


detectors shown in Figs. 4c, d,and e have single-lobe characteristics. The last detector can only be constructed with 


scintillation counters and belongs to the second type of compensated systems. It has a good geometrical symmetry 
since the main phosphor is located inside the compensating phosphor, 


On the basis of the above principles,we developed, in 1956-58, portable directional scintillation monitors for 


Y ~fay prospecting. * 


*Instruments using gas-discharge counters (type RGN-1) were also developed in 1956 (design due to V. Nilov). 
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The detector in one of the instruments (type RND-58) is in the form of a scintillation counter [FEU-31 and 
Nal(T1)], arranged in a compensated system of the first type. The main phosphor is placed in a lead filter 12 mm 
thick. The recording circuit is transistorized and is based on miniaturized components, The total weight of the de- 
vice is 3.5 kg, including a 1.5 kg probe. The instrument can be used to record y radiation in the range 10-104 pr/hr 
against a background of unwanted radiation having an intensity which is higher than that to be measured by a factor 


K 
a he 
se eo 
4K2 
Zz 


Zz z 
y Y y 
Ge x 
a) b) c) d) e) 


Fig. 4. Various designs of compensated radiation detectors and their 
directional properties. 


of 5. The directivity is 100-110°, and the introduction of an additional collimator reduces it by a factor of 5 to 8, 
The indicator is in the form of a differential ratemeter with a time constant of 3 to 4 sec and a microammeter. The 
"decompensation" does not exceed + (10-15)% when the spectral composition of the radiation changes by 100 kev. 
Total compensation can be obtained for any energy between 0.1 and 2 Mev. Satisfactory results have been obtained 
with this instrument under industrial conditions. 
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LETTERS TO THE EDITOR 


THE Li® (n, a)H? REACTION CROSS SECTION FOR 2.15 Mev 


NEUTRONS 


VoaP.. Perelygin and .K.. V - Tolstov 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 488-489, 
December, 1960 
Original article submitted March 3, 1960 


The total and differential cross section for the Li® (n, )H® reaction has been determined for neutrons from the 
D + D reaction having an average energy of 2.15 Mev. The reaction was studied in nuclear emulsion. 


Ilford-E-1 emulsions impregnated with lithium (containing 92% of Li®) were used. The lithium was introduced 
by soaking the plates in a 10% water solution of LiC,H,O, * 2H,O. The lithium concentration in the plates was ap- 
proximately (2.75 + 0.1) x 10" g/cm? and was determined by spectral analysis. A beam of singly ionized deuterium 
molecules having an energy of 500 kev was incident on a heavy ice target. The diameter of the target was 2.5 cm, 
The total neutron yield was determined with the aid of a proportional counter which recorded protons emitted in a 
known solid angle from the D(d, p)H® reaction. 


The emulsions were at a distance of 14.5 cm from the centerof the target,and the neutrons were incident nor- 
mally on the emulsion surface, Neutrons reaching the emulsions were at 145° to the deuterium ion beam, and their 
integral flux was (6.5 + 0.2) x 10’ neutron/cm”, The effective neutron energy was 2.15 + 0.1 Mev. At this energy 
the maximum range of the recoil protons in the emulsion was 45 ,while the minimum total track length of the @ 
particle and the triton was greater than 50. In addition to the 
recoil protons, protons having a range of the order of 654 were 
observed. These were due to the N“ (n, pCa reaction which 
takes place with the liberation of 0.6 Mev of energy. These 
protons constituted about 10% of the number of protons due to 
the reaction under investigation and could be easily identified 
by the lower ionization density and the absence of the a-par- 
ticle track, The ranges and angles of the a particles and the 
tritons to the direction of the neutron beam were determined in 
addition to their total energy. In some cases, owing to errors 
in the precise location of the reaction site in the nuclear emul- 
sion, the total energy differed by + 0.5 Mev from the true va- 
lue of 6.93 Mev. Only those events were selected for which 
there was a given relation between the range andthe angle. The 
number of cases for which this relation was not satisfied was 
less than 20% of the total number of measured tracks, 311 cases 
of Lif (n, &) H® reaction were observed over an area of 1.05cm’, 


46/a2, mb/sterad 


In the analysis of the results,corrections were introduced 
Differential cross section of the reaction Li(n, ce)H® for the yee’ leaving the emulsion, and the final calculation 
in the system with center of inertia Li’*. was foes out in the center of mass system of the compound 
nucleus Li, The reaction was not observed in the angle ranges 


0-15° and 150-180% owing to losses of tracks with large dip an- 
gles. 


. The differential cross section for the reaction is shown in the figure. The correction for the tracks having large 
uy angles was obtained by extrapolating the curve shown in the figure. The total reaction cross section, incorporat- 
ing all the corrections,was found to be 215 + 30 mb.* The root-mear-square error is made up of the following con- 


* This result was referred to earlier in ut 
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tributions: statistical — 6%, error in the number of nuclei — 4%, error in the determination of the neutron flux — 3%, 
geometrical corrections — 10%, spread in the values obtained in different series of experiments — 6%. 


Our result agreed to within experimental error with that reported in [2],where the reaction cross section was 
found to be 215 + 23 and 189 + 22 mb for neutron energies of 1.86 and 2.39 Mev. It is also in agreement with the 
results reported in [1], where the reaction cross section at 2.5 Mev was found to be 170 + 20 mb. The agreement is 
less satisfactory in the case of [3],where the reaction cross section for 1.5 and 2.0 Mev neutrons was found to be 320 + 
+ 60 and 270 + 40 mb,respectively. 


In conclusion, the authors wish to thank O, I. Kozints, T. A. Romanova,and F, A. Tikhomirov for help in the 
present experiments. 
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PRODUCTION OF FAST NEUTRINO BEAMS 


V. S, Barashenkov and Hsien Ting-ch'ang 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 489-490, 
December, 1960 
Original article submitted June 22, 1960 


The attention of many physicists has recently been turned to experiments with high-energy neutrinos [1]. These 
experiments give considerable possibility for investigations of many important characteristics of weak interactions. 
In particular, it is possible to check whether weak interactions actually become strong at high energies [2]. 


An intense beam of high-energy neutrinos can be obtained from large proton accelerators, such as the proton 
synchrotron at Dubna and the 30-Bev accelerator at Geneva. In this case, the neutrino is produced as the result of the 
decay of charged 7 mesons generated in collisions of fast nucleons (Fig. 1). The energy spectrum and angular distri- 
bution of such neutrinos can be calculated if the spectrum and angular distribution of the ™ mesons produced in NN 
collisions are known. In the general case, this is a very complicated problem, especially if the finite size of the neu- 
trino detector is taken into account [3]. However, the calculations are considerably simplified if only fast neutrinos 
emitted at small angles to the primary proton beam direction are considered. In this case, the angle of emission of 
the neutrino ©y differs very little from the angle at which the decaying 7 meson moves; Oy ~ © (see Fig. 1). 


Let us consider a neutrino of energy E > 1 Bev and angles Oy 
of the order of a few degrees. Since in the planned experiments 
the neutrino detector with a transverse cross section Sp ~ 1 m? is to 
be placed at distances of several tens of meters from the target, * 
such an approximation is quite satisfactory. 


If Wy (p; ©) is the momentum distribution of the charged 
mesons emitted at an angle © to the primary proton beam, the cor- 
responding energy spectrum of the neutrinos at a distance L from 
the target in which the 7 mesons are produced has the form 


W AED 6: Lia \ W,(p: 9) x 


>aE/m 


x [ 1—exp (- be =)] #R, : 
ap P 


Here, a = m? /(m?- 2), m™m is the ™-meson mass, and pt is the p-meson mass. If L is given in meters, then 7 =Tc- 107 2 
= 7.68 (T is the lifetime of a charged 7 meson). 


Fig. 1, Neutrino beam produced by the decay 
of ™ mesons generated by fast protons. 


This expression differs from the corresponding expression for y quanta produced in the decay of 7° mesons [4, 
5] only in the exponential term and the value of the coefficient a. 


Figure 2 shows the calculated neutrino spectra for the angle © = 0°. 


The calculations indicate that the spectra of neutrinos emitted at small angles depends very weakly on the 
angle, For example, the number of neutrinos of given energy emitted at an angle © = 3° is only 5% less than the 
number of neutrinos emitted at the same energy at an angle of © = 0°, This is also illustrated in Fig. 3, in which the 
values of the entire neutrino flux of energy E > 1 Bev and E > 2 Bev are shown for the angles 6 = 0° and © = 3°; 


ae! : 9 
Wy (8; L) = ( W,,(E; 8; L) dE; W°?) (@;L)= \ W,,(E; 8; L) dE. 
E>1 E>2 


“ 


We recall that, owing to the relativistic time dilatation, the ™ mesons decay only at large distances from the target. 
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The spectra Wy (E) and the functions wy) (L) shown in Figs. 2 and 3 are normalized to a single inelastic NN 
collision in the target and to a unit soild angle: 


\ W,, (E; @; L)dQdE=n, 


where n = 2, 3 is the mean number of charged ™ mesons produced in a single NN collision at an energy of 10 Bev 
[5, 6]. As in [5], the values of the function Wq (p; ©) was calculated from the statistical theory of multiple production. 


Fig. 2. Neutrino energy Fig. 3. Dependence of the neutrino flux of energy E> 
distribution at various > 1 Bev (1) and E > 2 Bev (2) on the distance L to the 
distances from the target target at neutrino angles of flight G = 0° (solid curve) 
L at an angle O = 0° to and O= 3°(dotted curve) with respect to the primary pro- 
the primary proton beam ton beam in the laboratory system. The fluxes Wy are 

in the laboratory system. normalized to a unit solid angle. 


The neutrino flux per unit solid angle Wy () (9; L) increases rapidly with the distance from the target. The 
limiting values of the flux for L= © are 


@ WS (8, co) W?) (@, co) 


0° 1,5 0,18 
a Ae 0,16 


With an increase in L, however, the solid angle of the neutrino 2p = Sp/ L? decreases rapidly; and,therefore, the num- 
ber of neutrinos recorded by the detector decreases with an increase in the distance between the detector and the tar- 
get. Thus, at a distance of L= 30 m,1.1 x 10~® neutrinos pass through 1 cm? of detector area (per single inelastic 
NN collision), and at distances of L equal to 50 and 100 m, the values are 0.6 x 10° and 0.27 x 1078, respectively. 
These values are very useful for various estimates in planning neutrino experiments. 


In conclusion, the authors express their deep appreciation to M. A. Markov for many discussions and valuable 
advice, and to R. Asanov and I. Polubarinov for discussions and valuable critical remarks, 
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A MODEL OF THE CIRCULAR SYNCHROC YCLOTRON 


V7 Aw Petukhov, I; "Gabanets, A-vA. -Churaviev. 

M. Karmasin, V. I. Kotov, E. A. Myaeé, Yu. L. Obukhov, 
Vo 2oknior, uw] sirakpre benda, I; Dobiashe Me Marek, 
Tt. Eukatko,sand. lL, Vess vetov 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 491-492, 
December, 1960 
Original article submitted May 28, 1960 


The idea of a circular synchrocyclotron, a strong-focusing accelerator with a time-independent magnetic field, 
was suggested in 1953 by A. A. Kolomenskii, V. A. Petukhov, and M. S. Rabinovich [1] and independently by Symon 
[2] in 1955, This type of accelerator would be very suitable from the point of view of increasing the intensity of the 
accelerated particle beam and the realization of the idea of colliding beams [3, 4, 5]. 


A model of the circular synchrocyclotron with radial sectors has been built at the Joint Institute for Nuclear 
Studies, The parameters of the model were chosen in such a way that it was possible to vary a number of its charac- 
teristics within a relatively wide range. In this way it was possible to carry out an exhaustive study of the behavior 
of the particle beam under various storage conditions. 


Figure 1 shows a general view of the accelerator. The electromagnet consists of eight periodic elements, each 
of which has normal and reverse sectors, and also of two straight-line sections. The azimuthal dimensions of the nor- 
mal sector, focusing the beam in the radial direction, are 22°, 30', while the dimensions of the reverse sector, which 
produces focusing in the vertical plane, and of the straight-line sections are 7° 30' each. The presence of the reverse 
sectors means that the orbit of the circular synchrocyclotron has a larger perimeter than that of an ordinary strong- 
focusing accelerator. This feacture of the circular synchrocyclotron is represented by a magnification coefficient 
which is equal to the ratio of the maximum radius of the. orbit to the minimum radius of curvature. In the present 
case this coefficient was apporximately equal to 3. 


The magnetic field in the sectors is produced by current coils wound on the surface of the poles along the radii 
of curvature of the instantaneous orbit. The coils are so distributed that the magnitude of the magnetic field increa- 
ses with the radius R of the orbit in accordance with the law H = Hg (R/Rp)*. The magnitude of the magnetic field is 
42 oe at R= 35 cm and 340 oe at R= 59 cm. In addition to the increase in the magnetic field with radius, the mag- 
net gap increases in proportion to R,so that the vertical dimensions of the working region vary from 2 cm at the initial 
radius to 4 cm at the final radius. The increase in all the geometrical dimensions of the sectors in proportion to the 
radius ensures the geometrical similarity of the instantaneous orbits. This condition, together with the requirement 
that the average field index k should remain constant (in the model k = 4), leads to the dynamic similarity of orbits 
[6, 7], i.e., to a constancy of the frequencies of free oscillations during the acceleration process. 


By changing the ratio of the fields in the normal and reverse sectors, it is possible to vary the number of beta- 
tron oscillations per revolution from 1 to 3 along the vertical, and from 2.5 to 3.5 along the radial direction. The 
fact that the frequencies of free oscillations could be varied in such a wide range means that it was possible to inves- 
tigate the effect of various resonances on the free oscillations, 


The model was designed for the acceleration of electrons. The injection of the electrons with energies between 
20 and 40 kev can be carried out continuously and also in pulses. The acceleration is produced by an induced elec- 
tric field (10-20 v per revolution) at a frequency of 450-500 cps. In order to improve the capture of the particles into 
the acceleration process, use is made of a forcing betatron core ensuring a voltage of 600 v per revolution over a time 
of about 1 psec. Experiments carried out on the model show that electrons can be captured into the accelerating pro- 
cess without the inclusion of the forcing core but the intensity in the absence of the latter is considerably lower. 
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Figure 2 shows a photograph of the accelerated beam with (a) and without (b) forcing. as this figure (c) ee 
a beam accelerated up to 450 kev, i.e., to the energy which can be communicated to the particles by the ed : 
electric field. By using a combination of the induced and high-frequency fields, the electron energy can be brought 
up to 2 Mev. 


Fig. 1. General view of the model of the circular synchrocyclotron; 


2) electron injector; 3) induction acceleration core; 4) high-frequenc 
multiplier, 


i) electromagnets; 
y oscillator; 5) photo- 


Preliminary results obtained at the initial stages of the operation of the model showed that the accuracy of the 


assembly and the stability of the main magnetic characteristics is extremely important. Moreover, in accordance 


with theory, a large number of various resonances were observed. These resonances have 


an important effect on the 
intensity of the accelerated beam, 
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Fig. 2. Pulses of the accelerated beam. a) 
Without forcing; b) with forcing; c) beam ac- 
celerated to final energy; d) magnetic flux in 
the main core. 
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APPLICATION OF THE SIMILARITY METHOD FOR SYSTEMATIZING 
EXPERIMENTAL DATA ON CRITICAL THERMAL FLUXES IN BOILING 


LIQUIDS 


S, S.*Kutateladze and G2 Il..Bobrovieh 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 493-494, 
December, 1960 
Original article submitted May 27, 1960 


Liquid-drop cooling is the most efficient method for cooling surfaces with intensive thermal fluxes. In this, 
the threshold of the maximum allowable thermal loads is determined by the value of the first critical thermal flux 
density for which film boiling occurs and the heat-transfer intensity sharply decreases (in any case, for a relative pres- 
sure p/ per Which is not very close to unity and at flow rates which are not too high). The large number of the para- 
meters influencing the value of the thermaFflux critical density makes the application of the similarity method highly 
desirable in this case. However, as is known, this method is effective only if the physical nature of the phenomenon 
under consideration is well understood [1]. 


At the present time, there are two hypotheses concerning the appearance of crisis under bubble-boiling condi- 
tions. One of them is most clearly formulated in papers by G. N. Kruzhilin [2 and 3]. According to this hypothesis, 
the first thermal-flux critical density is determined with respect to the same criteria which are used for the heat 
transfer processes in bubble boiling, namely: In the case of large-volume boiling in a liquid heated to the saturation 
temperature, the magnitude of the first thermal flux critical.density is determined by the criterion function [3] 


y" a [2 _cyAT" V 0 (y—y’) ee (1) 


eany 2, aay yale 


Te 
Vor AT! (y—Yy") x 


On the basis of many experimental data, represented in the form of products of the determining criteria in empirically 
chosen degrees, Kruzhilin proposed the following formula for calculations: 


= T'"0, 320,21 (y —y")0,48 (ry7)0138),0, 40 
Tor ~ 4700 ‘yor 31c0,08 40, 14 G (2) 


This hypothesis was not further developed. Nevertheless, Kruzhilin claims that it is valid from the physical point of 
view. 


The second hypothesis, which was put forward by one of the authors of the present communication, is based on 
the assumption that the crisis in boiling represents a qualitatively particular phenomenon which is connected with dis- 
turbances in the hydrodynamic stability of the two-phase boundary layer when the critical vapor generation rate is 
attained [4]. This hypothesis made it possible, by using the similarity method, to derive a number of equations which 
are accurate to a constant coefficient. Some of these coefficients were subsequently calculated by using theoretical 
methods [5]. 


For the case of large-volume boiling of a nonviscous saturated liquid, the hydrodynamic theory yields the equa- 


tion 


Ver 


ee =e COnS te 
r / os) Vo(y—y") = 


where K is of the order of 0,13-0.16. 


Thus, the hypothesis under consideration yields equations which greatly differ among themselves with respect to 
the characteristics of the boiling liquid physical properties figuring in these equations. At the same time, Eqs. (2) 
and (3) are in satisfactory agreement with a vast amount of data obtained in identical experiments. However, it should 
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be borne in mind here that Eq. (3) has been derived theoretically, while Eq. (2) is empirical in structure and has four 
empirical coefficients (the proportionality factor and the three exponents of the criteria under the function sign in 

Eq. (1)). Therefore, the agreement between experimental data and calculations based on Eq. (3) provides a more con- 
vincing confirmation that the process physical interpretation forming the basis of the derivation of Eq. (3) is correct. 


Complex k Values Calculated According to Eq. (4) BEER Sores reduced coitbe, fon 


aren ae cn 0.415 TM 3 (y — ymy0s28) 0-4 ; 
stance | ?/Pcr | k Istancel /Per ‘ yori y 20-14 7 0.64 ¢ 0,08, 0-04, 0.14 @ a r) 0.5 S 
Water! 0,005 0. 132 Ethyl ae . a we notice that it is a complex function of the saturation 
07086 | 0.121 |hol 07133 | 0,128 temperature and of a large number of the liquid physical pa- 
0,214 | 0,123 0,304 | 0,122 rameters. The table provides the values of the complex, 
ies men Hee wee which was calculated for a number of substances by using 
0,580 | 0,122 || Ben- 0,021 0,206 Eq. (4). It is obvious that the redundant quantities (more ac- 
Pro ete are zene en nae curately, the quantities which do not greatly affect the pro- 
pane 07390 | 0,187 0,390 | 0,180 cess under consideration), which appeared as a consequence 
ea as 0,670 | 0,214 of the incorrect. physical hypothesis, were combined in pro- 
0.86 0,227 cessing the experimental data into a complex which is in 
complete agreement with the hydrodynamic principles and 


which remains a practically constant quantity. 


Thus, Kruzhilin’s unsuccessful approach in attempting to derive an equation for the critical thermal flux leads, 
even in the simplest case, to the introduction of a large number of factors which only slightly affect the process in 
question and, correspondingly, to an unnecessary complication of the equation. This was also the reason why Kruzhilin's 
method did not bring about any further advances in studying the effect of forced motion, underheating of the flow 
center, which makes it impossible to attain the saturation temperature, and other factors on the critical phenomena 
in boiling liquids. 
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HEAT TRANSFER TO A MELT OF SODIUM AND POTASSIUM 


IN ANNULAR CLEARANCES 


E. M. Khabakhpasheva and Yu. M. Il'in 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 494-496, 
December, 1960 
Original article submitted January 14, 1960 


The majority of papers on heat transfer to liquid metals which were published in recent years in our country 
and abroad were devoted to the investigation of heat transfer in round pipes. Heat transfer in channels with more 
complicated sections has been studied to a much lesser extent. The theoretical solutions obtained in determining the 
heat-transfer coefficients in the case of flat and annular clearances [1 and 2] are based on a number of simplifying as- 
sumptions, and they must be checked experimentally. 


In most cases, experiments on heat transfer in annular clearances were performed in different heat exchangers, 
and the values of the heat-transfer coefficients were determined by calculation on the basis of the measured heat-trans- 
fer coefficients [3]. This work involved the use of various assumptions concerning the relationship between the heat- 
transfer coefficients at the inside and the outside surfaces of the dividing pipe. 


In the present paper, the heat-transfer coefficients in an annular clearance were determined on the basis of tem- 
perature measurements in a section sufficiently removed from the channel inlet. The experimental section consisted 
of a thick-walled copper tube with an inside diameter of 17 mm, Heaters (or flow restrictors) were installed along the 
pipe axis and were centered by means of special bushings at the experimental section inlet and outlet. The wall tem- 
perature along the copper pipe length was measured by means of eight thermocouples. In order to eliminate the heat 
leakage from the heater, a thermometer and a compensating heater were installed in the experimental section. For 
controlling the liquid-metal purity and the heat-transfer surface cleanliness, the heat-transfer coefficient was periodi- 
cally measured in a pipe with a round cross section (without the internal heater). .The heat-transfer coefficient values 
in the pipe were stable and were in good agreement with Lyon's equation [4] (Fig. 1): 


Nu=7-+0,025 Pe®8, 


The experiments on heat transfer to a melt of sodium and potassium were performed for two-sided and single- 
sided heat supply in annular clearances of different widths (2.5, 3.5, and 4.5 mm), The experimental results are shown 
in the form of graphs in Figs. 2 and 3, where the Nusselt and Peclet numbers were obtained on the basis of the equiva- 
lent hydraulic diameter de = dy-d2, which was used as the characteristic dimension. The diagrams also show curves 
representing the results of theoretical calculations [1 and 2]. 


It follows from Fig. 2 that, for one-sided heat supply, the Nusselt number values for narrow annular clearances 
(d\/d2 © 1) are lower by 25 to 30% than the Nusselt number values for round pipes for the same Peclet number values, 
and that they are in satisfactory agreement with theoretical data. The Nusselt number values for annular clearances 
and for two-sided heat supply (see Fig. 3) in cases where the Peclet number exceeds 500 also are in satisfactory agree- 


ment with theoretical data [1 and 5]; and they are 1.3 to 1.5 times as large as the corresponding Nusselt number va- 
lues for round pipes. 


Thus, the above results confirm the theoretical conclusion that, for coolants with small Prandtl number values 
changes in the channel geometric shape or in the method of supplying heat considerably affect the heat-exchange 
intensity, The recommendation contained in manual [3], where it is suggested that the heat-transfer coefficients for 
liquid metal flow in flat and annular clearances be calculated in the first approximation by using heat-transfer equa- 


tions for pipes with round cross sections on the basis of the equivalent hydraulic diameter,can be considered only as 
a very coarse approximation. 


, It should also be noted that the theoretical recommendations can be applied in heat-transfer calculations only 
if no thermal contact resistance is present. This condition is apparently satisfied in the case of heat transfer to alkali 
metals if they are thoroughly freed from oxides [6]. It is presently not yet possible to eliminate the thermal contact 
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300 400 500 600 700 «= 800——« 900s «000s tt00—s «Pe 
Fig. 2. Heat transfer to Na — K in annular clearances with different widths for one- 
sided heat supply (d, is the diameter of the heated surface). Nupj = 7+ 0.025 Pe’® [4] 
is the Nusselt number value for round pipes. 


5+0,0360°8[ 5) 
® aay 
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Fig. 3, Heat transfer to Na — K in an annular clearance for two-sided heat supply 
[(dy/ da) = 1.4; (dy = d2)]. Here, qy and qp are the outside and inside thermal fluxes. 
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resistance in the case of heavy metal flow (mercury and the eutectic lead and bismuth alloy) in steel pipes [7]. In 
these cases, the effective transfer coefficient values for pipes with round cross sections and annular clearances depend 


on the thermal contact resistance value. 
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THERMAL CONTACT RESISTANCE 
Yul Po shiykov and E. A. sGanin 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 496-498, 
December, 1960 
Original article submitted May 28, 1960 


The thermal contact resistance must be taken into account in the case of large thermal flux values [1 and 2]. 


The actual contact area constitutes only an insignificant portion of the over-all area of the surfaces in contact 
(in most cases, it does not exceed 1-2%[3 and 4]). The heat transfer between surfaces in contact occurs mainly due 
to heat conduction through the medium filling the space between the burrs in the contact zone and through the points 
of direct contact. 


The over-all thermal resistance can be expressed by the relation 


ee Ee 
Roa Rgas HKact ’ (1) 


where Rog is the over-all thermal contact resistance, Rgas is the gaseous layer thermal resistance, and Ract is the di- 
rect contact thermal resistance. 


For Rgas, we can use the thermal resistance of a layer whose thickness is equal to the equivalent (with respect 
to the volume) clearance between the surfaces in contact: 


Vg 
Sp (2) 
where Vegas is the volume enclosed between the projections in the contact zone, which is filled with a gas whose ther- 
mal conductivity is Agas, and Sp is the contact nominal area. 


Since the absolute value of deg is negligible for the classes of finish quality accepted in machine construction, 
the heat can be transferred through the layer only by pure heat conduction, i.e., 


gas A ‘ (3) 
Let us use the concept of the relative drawing together of surfaces in contact when a load is applied [5]: 


2 iad (4) 


where a is the absolute value of the reduction in spacing (1), and hay is the microroughness height (1). 


By expressing Vgas and then Seg in terms of e for the case where a rough and an “absolutely” smooth surface 
are brought into contact, the following relationship is obtained for Rgas: 


h 
Se 
Ryas = at (1—e)?. (5) 


For determining Ract, we used the results obtained in an analytical investigation of contact spot thermal resis- 
tance, which was performed by Cetinkale and Fishenden [6], and the results obtained in Holm's work on electrical 
contact resistance [7]. 


When surfaces are compressed, the contact spot dimensions remain practically constant, and the increase in the 
actual contact surface area occurs mainly as a result of the greater number of contact points [1, 5, and 8]. If we de- 
note by 7 the ratio of the actual contact area Sact to the nominal area Sn, we obtain the following dependence: 
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Risk te ed (6) 
act  21ykm 3 
where Am is the thermal conductivity of the material in contact, and 2.1 x 10° is a coefficient whose dimension is 
1/m. 

The relative magnitude of the actual contact area 7 and the relative reduction in spacing € vary in relation 
to the load and the quality of the surface finish. These factors determine the thermal resistance variation. 


On the basis of theoretical and experimental investigation of contacts between rough and absolutely smooth 
surfaces, the relations of ¢ and 7 to the load, the microroughness height (finish quality), and the physical prop- 
erties of materials, were obtained in [5]. 


By using expressions (5) and (6), the value of Roa can be de- 
termined with respect to known values of € and 7. 


Calculations show that the value of (1-e) and, consequently, 
of Rgas does not change by more than 15% (Fig. 1); and, therefore, 
for practical purposes, Rgas can be determined by using the simpli- 
fied equation 


Royas = 97 (7) 

For surfaces whose finish quality is not better than class 9 or 
10, the following calculation formula was obtained after some simp- 
lifications [9]: 


SpCo. 


os 24-108 kN” (8) 


0,2 4 
Oa exper. < 
es where N is the load (kg), Os is the yield point of the material after 
it has been subjected to the maximum amount of cold-working (it 
is assumed to be equal to the ultimate strength), and C is a coeffi- 
cient which is assumed to be equal to 3. 


50 100 150 P, kg/cm? Figure 1 shows the theoretical and experimental curves for the 


Fig. 1. Dependence of the thermal resistance = case where a smooth (V 12) surface is in contact with a rough (V4) 
and its components on loading. surface. The specimens were made of Steel 1Kh18N9T. 


An analysis of the Rgas theor and Ract theor curves indicates that the over-all thermal resistance components 
are commensurable within a rather wide loading range. The decrease in the over-all thermal contact resistance when 
the specimens are compressed with a greater force is mainly due to an increase in the actual contact surface area. 


The contact between two rough surfaces is usually encountered in practice. In this case, Rgas is calculated by 
using the equation 


d gas (9) 
where Ract is determined from expression (8). 


The results obtained in analytical calculations of thermal contact resistance according to Eqs. (8) and (9)* are 
in good agreement with experimental data. Some of these data are shown in Fig. 2. 


Dinah 


hin MmitAme * 


7 


Persie ; 
If different metals are in contact, the 0, value for the less plastic material is used, and Ay, is determined from the 
expression shown above. 
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Thus, the thermal resistance of the contact between rough surfaces for different compressive forces can be de- 


termined analytically with a sufficient degree of accuracy. 


PAP LY 


fe (m2. hr -°C)/kcal 


g 50 100-150 kg/cm? 


Fig. 2. Thermal contact resistance. 1) Steel 
1Kh18N9T, finish; V5; 2) duralumin D-16, 
finish: V4; ----) theoretical data. 
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EFFECT OF RADIOACTIVE RADIATION ON THE DIELECTRIC 


PROPERTIES OF ELECTRICAL INSULATORS 


K. A. Vodop'yanov,-B. I. Borozhtsov, M. D. Lavrov, 
E. S. Nesmelova, and G. I. Potakhova 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 498-500, 
December, 1960 
Original article submitted December 1, 1959 


The investigation of the effect of y radiation on the dielectric properties of solid plastic dielectrics is of great 
importance for the selection of electrical insulators for electrical and electronic devices which are exposed to radia- 


tion. 


In dependence on the physicochemical properties of the material and its structure, radiation can cause ioniza- 
tion, displacement of atoms, or dissociation of molecules. In this, the electrophysical properties of solid dielectric 
changes considerably. 


We investigated the frequency and the temperature characteristics of the permittivity and the dielectric loss 
angle for polyethylene, “Kel-F," and “product 10" (a mixture of polystyrene with vinylnaphthalene) before and after 
they were irradiated by y rays; the dose rate was from 400 to 1200 r/min, and the irradiation dose was from 2000 to 
100,000 r, A betatron with 15-Mev particles served as the y-radiation source. The specimens consisted of disks which 
were 1-2 mm thick. 


The frequency dependences of the loss angle and the per- 
mittivity of the organic dielectrics under investigation remain 
almost unchanged after irradiation with a dose of 50,000 r. For 
a frequency of 10" cps and an irradiation dose in excess of 
1000, 000 r, the loss angle slightly increases only in "Kel-F"; 
for polyethylene and the “product 10,” the loss angle variation 
remains within the measurement error limits. The fact that 
the loss angle is independent of the radiation effects in all of 


20 40 CO 80 100 f20 140 160 180 200 


7% the three dielectrics is also observed at different temperatures 
Fig. 1. Dependence of the loss angle on tempera- DE 
ture at the 10 cps frequency for an irradiation dose The permittivity of polyethylene, "Kel-F,” and “product 
of 50,000 r, 1) Polyethylene; 2) "Kel-F"; 3) “pro- 10" changes little after they are exposed to temperature action 
duct 10"; O) unirradiated; x) irradiated. and radiation, On the basis of an analysis of the frequency and 


the temperature dependences of the loss angle and permittivity 
of organic dielectrics before and after irradiation, it can be assumed that, in polyethylene, the frequency and the tem- 
perature maxima of the loss angle are determined by directional polarization of the polar radicals (c = 0) in molecules, 
With respect to its properties and molecular structure, “Kel-F" is close to polyethylene, with the difference, however, 
that "Kel-F" contains fluorine, which is characterized by stable bonds, instead of hydrogen with unstable bonds. “Kel- 
o is a weakly polar dielectric, and, therefore, the observed increase in the loss angle with an increase in frequency 
is apparently also due to the orientation of polar groups entering the "Kel-F" composition. In comparison with poly- 
ethylene and "Kel-F," the loss angle in “product 10" depends to a much lesser degree on frequency and temperature 
which indicates that dielectric dipole losses are not present in the measured frequency and temperature intervals. 


Figure 2 shows the dependences of the loss angle and the permittivity of glass textolite SKM-1 on frequency 
before and after irradiation. In order to determine the causes of the increase in the loss angle of glass textolite 
SKM-1 after irradiation in the low-frequency region, we measured the temperature dependence of the specific vo- 
lume resistance (Fig. 3) and the loss angle of glass textolite SKM-1 at different frequencies (Fig. 4). On the basis of 
the obtained data, it can be assumed that the effect of radioactive irradiation on glass textolite SKM-1 is manifested 
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Similar characteristics were obtained for the AG 4, K-211-3, K-114-35, and FKIM-25 plastics, which are based 
on various modifications of phenolformaldehyde resins. In dependence on external conditions (temperature and fre- 
quency), the dielectric loss angle in these plastics will be determined by relaxation processes and by conductivity 
(Fig. 5). Loss angle measurements, performed at a temperature of -60°C, indicate the absence of conductivity losses 


and the presence of relaxation losses. 


item es 4 Si 6 87 8 oof 


Fig. 2. Dependence of the loss angle and the 
permittivity of glass textolite SKM-1[(O) un- 
irradiated; x) irradiated)]on frequency atroom 
temperature for an irradiation dose of 100,000r. 
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Fig. 3. Dependence of volume resistivity 
of glass textolite SKM-1[(O) unirradiated; 

X) irradiated)] on temperature for an irra- 
diation dose of 100,000 r. 


The temperature dependences of the loss angle in polyamide- 
68 at different frequencies are shown in Fig. 6. These dependen- 
ces can be used as an example for explaining the manner in which 
irradiation affects the various losses. The decrease in the angle 
loss, which is determined by and is observed at frequencies from 
50 to 10 cps, is probably connected with an increase in the poly- 
merization of polyamide-68 under the action of radiation. The 
fact that dipole losses did not change after irradiation is connected 
with the prevalence of the orientation effect of polar groups in 
polyamide-68 over those changes in the dynamics of elementary 
particles or molecules which are caused by irradiation. Similar 
phenomena are observed in other polar plastic dielectrics, such as 
polyvinyl chloride, lavsan, and FKPM-25, 


We also investigated two varieties of organic silicon rubber: 
rubber 14r-2, which is based on synthetic organic silicon rubber, 
titanium dioxide, and benzoyl peroxide, and rubber 141-6, which 
differed from rubber 14r-2 by the fact that it contained powdered 
silica gel and zinc white instead of titanium dioxide. In compari- 
son with rubber 14r-6, rubber 14r-2 had a lower electrical conduc- 
tivity and a smaller dielectric loss angle. 


The results of loss angle vs frequency measurements before 
and after irradiation indicate (Fig. 7) that a decrease in the loss 
angle with an increase in frequency is characteristic of both rub- 
bers, The effect of radiation on rubber 14r-6 depends on the state 
of the specimen under investigation. If the specimen has not pas- 
sed through the second vulcanization stage, as was the case in our 
experiments, irradiation leads to a sharp decrease in the loss angle, 
which diminishes with an increase in frequency. It is obvious that 
irradiation leads to further vulcanization of rubber 141-6. A dif- 
ferent situation is encountered in the case of completely vulcanized 
specimens; After such specimens are irradiated, the dielectric loss 
angle and permittivity increase, while the effect of irradiation is 
similar to that observed in rubber 14¥-2, 


The presence of ohmic dielectric losses in rubber apparent- 
ly also causes a noticeable change in the dielectric loss angle in 
the low-frequency range after irradiation. 


Thus, for the majority of solid plastic dielectrics, the change in the loss angle as a result of irradiation is the 
greater, the lower the frequency. With an increase in frequency, this change diminishes and becomes hardly notice- 


able at high frequencies. 


We are of the opinion that the loss angle variation in solid plastic dielectrics as a result of radioactive irradia- 
tion is determined by the dielectric loss mechanism, It was shown experimentally that, in cases where ohmic losses 
are observed, also a considerable change in the loss angle occurs; in the presence of dielectric dipole losses, notice- 


able changes in the loss angle do not take place. 
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Fig. 4. Dependence of the 
loss angle in glass textolite 
SKM-1[(O) unirradiated; x ) 
irradiated)] on temperature 
for an irradiation dose of 
100,000 r at the following fre- 
quencies: 1)106 cps; 2 and 
210, cps; 3 and 3")50 cps. 


Fig. 5. Dependence of the loss angle in plastic AG-4 
((O) unirradiated; x) irradiated) on frequency for an 
irradiation dose of 48,000 r at the following tempera- 
tures: 1) + 20°C; 2) -60°C. 
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Fig. 6. Dependence of the loss angle in polyamide- 
68 [(O) unirradiated; x) irradiated)] on temperature 
for an irradiation dose of 50,000 r at the following 
frequencies: 1 and 1") 50 cps; 2) 106 cps; 3 and 3") 
5 + 10° cps. 


Fig. 7. Dependence of the loss angle in rubber [(O) un- 
irradiated; x) irradiated)] on frequency for an irradia- 
tion dose of 50,000 r. 


Consequently, before any conclusions are drawn re- 
garding the resistance of solid plastic dielectrics to the ef- 
fect of a certain dose of y -radiation of a given intensity, 
it is necessary to know the character of polarization and 


of dielectric losses in different temperature and frequency 
intervals. 


The authors hereby express their gratitude to N. I. 
Ol'shanskaya, T. G. Mikhailova, L. T. Murashko, and A.1. 
Tovbina for their help in performing the experiments. 


ISOTOPIC COMPOSITION OF RUTHENIUM 
K. G. Ordzhonikdze and O. S. Akirtava 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 501-503, 
December, 1960 
Original article submitted May 4, 1960 


The relative content of ruthenium isotopes has been determined on several occasions [1-3]; but the method 
employed did not have the simplicity which is desirable for series analyses, and in a number of cases did not make it 
possible to carry out precision measurements. The problem was therefore presented of developing a fairly accurate 
and simple method of measuring the isotopic composition of ruthenium, which does not require preliminary chemical 
treatment of the sample and permits direct measurement of the isotopic composition of enriched samples, taken from 
the pockets of an electromagnetic separating plant. 


The experiments were carried out in an MS-3 type industrial mass spectrometer with recording of the ion cur- 
rents by an external microammeter and a Speedomax Recorder automatic electronic potentiometer with magnetic- 
field linear time scanning. A standard source with focusing of the ion beam was used. Check measurements in neon 
and mercury (Hg* and Hg**) showed that systematic errors do not occur if the apparatus is correctly adjusted [4, 5]. 


To obtain maximum accuracy of the measurement of the isotopic composition of ruthenium, the conditions for 
obtaining the most stable ion currents during the ionization of ruthenium vapor by an electron beam and by the ther- 
mionic emission method were investigated. In the first case, a directly heated evaporator, making it possible to 
avoid overheating of the components of the source, was used to obtain ruthenium vapor. The inverted U-shaped evap- 
orator, made of 2 x 0.02 mm tungsten strip, was constructed to a width of 0.9 mm in the center and had a cavity, 
holding 1-2 mg of the substance investigated. This shape made it possible 
to obtain the temperatures necessary for the evaporation of ruthenium, with 
the minimum heating power (~ 30 watts). The evaporator was fixed in the 
holder installed opposite ‘the corresponding slit in the ionization chamber 
of the ion source. The design of the holder makes it possible to change the 
evaporator rapidly. To increase the stability of evaporation, ruthenium 
powder was mixed with degassed tungsten powder in the optimum ratio of 
1:3. When this mixture is sintered it becomes porous and evaporation of 
ruthenium from it takes place more uniformly; the process continues longer 
and the ruthenium evaporates completely, without leaving a residue. With 
ruthenium ion currents sufficient for measurements, this made it possible to 
reduce the weight of the investigated sample for analysis. Thus, to carry 
Fig. 1, Emitter with holder; 1) Emit- out an experiment 0.3-0.7 mg of ruthenium metal was placed in the evap- 


ter; 2) ceraniic insulator of the holder, Orator. 


Depending on the character of the source, two maxima of the Ru‘ ion current were observed with corresponding 
electron ionization energy value of 20 and 90-100 ev, respectively. In the latter case the line strength of the ruthe- 
nium isotopes was two-threefold more intense than with an electron energy of 20 ev. In addition to the monovalent 
ions, di- and trivalent ruthenium ions were observed in the spectrum of the masses. The value of the current of Rutt 
ions was three- fivefold less intense than the current of Rut ions, and the emission intensity of the Ru**t ions did not 
allow them to be measured in practice. Lines, corresponding to ruthenium compounds, were absent in the observed 
spectrum, 


Measurements of the relative abundance of ruthenium isotopes were carried out in currents of univalent and di- 
valent ions. To reduce the mercury background, a series of experiments with an electron ionization energy of 20 ev 
was carried out. The current of Hgtt ions was reduced 20-30-fold, so that the background was reduced to 0.1-0.5%, 
of the useful signal, Within the limits of such accuracy the background could be neglected. 


Measurements in divalent ruthenium atoms were complicated by the presence of background lines with mass 
numbers of 50, 51,and 52 in the spectrum. The background was due to the material of the evaporator, and its inten- 
sity reached 5% of the value of the useful signal. 
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By means of the method of ionization by evaporation and subsequent bombardment by an electron beam, the 
relative content of ruthenium isotopes can be determined with an accuracy reaching 0.5-1%, Within the limits of 
such an error,it is best to carry out the measurements in Rut ions with a minimum energy of the ionizing electrons 


(20 ev), when the mercury background becomes negligible. 


To obtain a high accuracy of measurement of the isotopic composition of ruthenium, we investigated we passe 
bility of its thermionic emission; for this purpose, an emitter (Fig. 1),which was previously employed for the isotopic 
analysis of alkali metals [4, 7], was used. This emitter was punched from tungsten foil of thickness 0.02 mm. The 
area of its central part (14 x 0.9 mm) corresponds to the dimensions of the slit of the first lens of the ion source. The 
U-shaped profile of the active part does not permit appreciable deformations of the emitter, which could influence 


the form of the lines of the mass-spectrum. 


A mixture of ruthenium powder and degassed tungsten in a gravimetric 
ratio of 1:5 to 1:20 was wetted with alcohol, and about 6-8 mg was placed in 
the cavity of the emitter, forming an even layer. 


; i iy ste -14 
An appreciable emission of ruthenium ions, of the order of 10 


for Ru!” isotope, commenced at an emitter temperature of ~ 1950°C. 


amp 


Measurements of the abundance of ruthenium isotopes were carried out 
with currents (at the collector) of the order of 10712-10714 amp in the 2000- 
2250°C temperature range. The pressure on the source did not exceed 5 - 107" 
mm Hg, and remained at(1-2) - 107? mm Hgon the analyzerside, Variations of 
intensity within the limits of sensitivity of the apparatus were not observed. 
The value of the current of ruthenium ions remained constant for a certain 
time after the commencement of emission, depending on the amount of 
sample used, and then decreased regularly as the material was exhausted. The 
intensity could be kept constant by gradually increasing the heater current. 


=? a 


Figure 2 gives the mass-spectrum of ruthenium isotopes obtained as a 
result of the observations (the sensitivity of the apparatus was increased 100% 
56 58 39 100 101 102 104 for recording the Rut and “Rut peaks). To obtain the same measurement 
accuracy for all isotopes, the sensitivity of the recording circuit was increased 
during the recording of rare isotopes, so that the deviation of the potentio- 
meter carriage was the same for all isotopes. The error of the determination 
of the relative abundance of each ruthenium isotope did not exceed 0.2%, 
To exclude the influence of the diminution of the intensity on the measurement results, scanning was carried out al- 
ternately in the direction of an increase in the masses and in the opposite direction. The abundance of the isotope 
was determined for one experiment as the mean of an even number of recorded spectra (in our case, 10). Table 1 
gives the measurement results obtained in this way. 
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Fig. 2, Mass-spectrum of ruthenium 
isotopes. 


As is shown by the data of Table 1, the statistical errors in the determination of all isotopes, with the exception 
‘of *Ru, are within the limits of 0,06-0.15%. The value of the accuracy of the determination of the relative content 
of each isotope, given by the recording circuit, i.e., 0.2%, is given as the lower limit of errors for the values found 
for the relative abundance of ruthenium isotopes. **Ru isotope, for which the statistical error somewhat exceeds this 
value, forms an exception, 


The values we obtained for the relative content of ruthenium isotopes are given, together with the results of 
other authors, in Table 2, As may be seen from this Table, the results of our experiment agree closely with the va- 


lues given in [3]. The data of [1-2] are somewhat too high for the rare isotopes and correspondingly low for the more 
common isotopes. 


In spite of the high ionization potential of ruthenium (7.3 ev), considerably exceeding the work function of 
pure tungsten (4.52 ev), a fairly high value is obtained for the emission intensity of the ruthenium ions (Lontasort 
at the collector of the apparatus), which is also stable, ensuring high accuracy (0.2%) of isotopic analysis. In all prob- 
ability this can be explained by the fact that the work function of some parts of the sintered tungsten powder was 
increased as a result of the monoatomic layer of oxygen adsorbed on them, But the surface film of oxygen is very 
stable up to temperatures of 1700-2200°C [6]. The ionization of the ruthenium atoms in these parts of the emitter 
evidently also ensures a sufficient emission intensity for accurate measurements. 
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TABLE. 


Data of the Mass-Spectrometric Determination of the Content of Ruthenium Isotopes 


| 


Expt. No, 96Ru 9sRu 9Ru | 100Ru | 101Ru | 102Ru | 104Ru 

1 5,432 1,869 12,62 12,52 16,95 31,97 18,63 

2 5,418 1,885 12,95 12,46 17,04 31,74 18,92 

3 5,442 1,883 12,66 472753 16,99 31,63 18 ,90 

4 5,392 1,861 42,62 12,54 17,03 31,68 18 ,88 

5 5,494 1,884 12,64 12,60 17,14 31,37 18,87 

6 0,476 1,884 12,60 12,53 17,05 34,61 18,85 

qu 53,490 1.,864 12,61 12,54 17,06 34,68 18,76 

8 5,479 1,858 12,66 12,54 17,04 34,64 18,82 

9 9,471 4,856 12,62 12,54 ASE) 34,66 18,86 

10 3,409 1,876 12,65 12,94 16,98 31,93 18,98 

14 5,486 1,867 42,65 12,57 ale (OK 31,59 18 83 

12 5,900 1,880 12,64 12,53 17,03 31,48 18,94 

13 9,472 1,816 12,62 12,47 16,96 31,64 19,02 

Mean 5, 460 1,868 12,63 12,53 17,02 31,63 18 ,87 
Mean | 

square 
error | 0,008 0,005 0,008 0,010 0,014 0,039 0,027 
TABLE 2 


Results of Determination of the Isotopic Composition of Ruthenium from the 
Data of Various Authors 


Values obtained 


Ruthenium | | | ; 

j [1] [2] [3] in the present 

isotopes work 
*®Ru 5,68+0,17 | 5,50+0,02 | 5,47+0,03 | 5,46-+0, 04 
Ru 2,22+0,09 | 1,9140,01 | 1,840,02 1,868-£0, 005 
RU 12,810,09 | 12,700,041 | 12,770, 05 12,63-£0,02 
200 Ru 12,70+0,19 | 12,69+0,05 | 12,56--0,05 12,53--0, 02 
101Ru 16,98+0,25 | 17,01+0,03 | 17,10--0,06 17,02--0,03 
102Ru 31,34+0,47 | 31,52+0,06 | 31,70-£0,06 31,6340, 06 
104Ru 18,270.27 | 18,67-£0,06 | 18,56-L0,05 18 ,87-F0,04 


In all probability, the use of a source with separate evaporator and ionizing filament [8] for ruthenium metal 


could not give an appreciable increase in the ionizing efficiency because the temperatures at which thermionic emis- 
sion of ruthenium is observed are not markedly different from the maximum working temperature of the ionizing fi- 
lament. Moreover, a sufficient vapor pressure of ruthenium metal is obtained only at high temperatures, close to the 


ionization temperature. 
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THE HEAT OF FORMATION OF PuBejs 


V. V. Akhachinskii and L. M. Kopytin 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 504-505, 
December, 1960 
Original article submitted July 5, 1960 


The literature contains no data on the heats of formation of intermetallic plutonium compounds. 


We have determined the heat of formation of PuBej3 by measuring the heat of solution of this intermetallide and 
its components in 19% hydrochloric acid, using an isothermally jacketed microcalorimeter (see figure), The calorim- 
eter vessel consisted of an inner container welded from sheet tantalum and an outer container turned from copper. A 
manganin wire heater of 180 ohms resistance was placed in the space between these containers and the region filled 
with a layer of paraffin to hermetically seal the inner container. A copper resistance thermometer (40 m of wire, 
0.05 mm in diameter) wound on the outer surface of the calorimeter was incorporated in a bridge circuit containing 
a lamp and scale galvanometer, the latter having a sensitivity of ~ 0.00003° per mm scale deflection. The thermal 
equivalent of the calorimeter was approximately 35 cal/deg and the cooling constant, 3 - 10-* sec”*. The tempera- 
ture of the jacket water (25°C) was held constant to within + 0.001° by a special electrical device [1] which utilized 
a MMT-9 thermistor of ~ 1000 ohm resistance as a sensor. Distillation of the solvent was avoided by so choosing the 
temperature at the beginning of the principal period that the temperature at the final period would be somewhat be- 
low the jacket temperature. The calorimeter stirrer was lowered several millimeters at the beginning of the principal 
period to break the container on a special protuberance on the bottom of the calorimeter and bring the working ma- 
terial into reaction with the solvent. The hydrogen evolved during the experiment was collected in a gas burette. 
The calorimeter was cooled to 22° at the end of dissolution and its thermal equivalent determined by introducing a 
definite quantity of electrical energy. The results are expressed in 15° calories. 


A heat-exchange correction was obtained from the formula: 


n—t1 
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sai( 5, 


Joe ene, _ i=1 4 


i 


in which k represents the cooling constant for the calorimeter, Sc is the convergent temperature, 9p is the tempera- 
ture at the end of the initial period, $, is the temperature at the end of the principal period, and At, is the time in- 
terval between the last reading in the initial period and the first reading in the principal period. 


Measurement of the amount of liberated hydrogen (with the requisite corrections) made it possible to determine 
the chemical composition and the phase composition for each charge of alloy from the well-known phase diagram for 
the Pu-Be system [2]. The phase composition was calculated from the equations 

+ y= 1; 
138 802 +-908,26y =v, 
in which x and y are the amounts of Pu and PuBe4s, feppectively. per 1 g of alloy; the coefficients 138.80 and 908.26 
represent ‘the amounts of hydrogen expressed in STP cm? evolved in the dissolution of 1 g of Pu and 1 g of PuBejs, res- 
pectively (the latter figure was obtained from data on the component Pu and Be); and v is the volume of hydrogen 
evolved in the dissolution of 1 g of alloy. 


The Be used in this work was 99.15% pure. The plutonium contained less than 0.2% total contaminants. The 
intermetallide PuBe,3; was prepared by fusing a mixture of metallic plutonium chunks and powdered beryllium in a 
high frequency electrical furnace. Fusion was carried out in a BeO crucible in an atmosphere of nitrogen. The lat- 
tice constant of the product compound (face-centered cubic) was a = 10. 259 + 0.001 kX* and its microhardness, 


* The lattice constant of PuBeys, a, is 10.26 kX according to [4]. 
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Results on the Dissolutions of Beryllium, Plutonium, and the Alloy 


Dissolved 
substances 


Number of 
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Construction of microcalorimeter: 


1) Calorimeter 


jacket; 2) thermistor; 3) Beckmann thermometer; 4) 


housing; 5) calorimeter vessel; 6) container with working 


material; 7) impellers; 8) calorimeter cover; 9) ebony 


bushing; 10) tube; 11) felt packing; 12) copper. bushings; 


13) tube for escape of hydrogen; 14) packing nut; 15) 


spring; 16) Warren motor; 17) stirrer shaft; 18) tube con- 


necting calorimeter with gas burette; 19) jacket heater; 


2) jacket cover; 21) stirrer. 
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H, evolved per 1g of 


experiments substance, c 
ured at STP) 


Heat of so- 
lution, AH, 
kcal/ mole 


m°( meas 


89.38 + 0.06 
141.022 0.19* 
1267.24 2.3 


9477.8 + 0.67 
138.8 + 0.13 
From 830 to 848 


2 


1045 kg/mm*. Metallographic analysis showed small in- 
clusions of free plutonium distributed on the grain bounda- 
ries of the intermetallide. Measurements on the amount of 
hydrogen evolved during alloy dissolution showed this plu- 
tonium to be distributed nonuniformly. The phase compo- 
sition of the charge was determined from these same meas- 
urements and varied from 92.17 weight 7% PuBe,3 and 7.83 
weight % Pu to 89.78 weight % PuBe,3 and 10.22 weight 

Yo Pu. 


The heat of solution of PuBe,3 was calculated for 
each charge separately, account being taken of the phase 
composition in order to avoid effects from nonuniformity 
of the alloy. 


Results on the dissolution of beryllium, plutonium, 
and the alloy are presented in the table. The values for the 
heat of solution contain a correction for the heat of vapo- 
rization of the hydrochloric acid required for saturating the 
liberated hydrogen (0.206 kcal/mole H ). 


The probable error is given as a measure of the ac- 
curacy of the results for Be and Pu. The maximum expe- 
rimental error is given for PuBes3, this quantity covering 
both accidental errors and possible systematic errors in the 
determination of the alloy composition. These latter er- 
rors arise from the fact that the quantity of gas evolved in 
the dissolution of unit mass of a mixture of the components 
is not the same as that evolved in dissolving the alloy ob- 
tained from this mixture. This, in turn, is due to such fac- 
tors as the vaporization of the volatile contaminants and in- 
teraction of the contaminants with one another. 


The data of the table were used for a Hess Law cal- 
culation of AH}9g, the heat of formation of PuBe3, and 
led to the value of 35.7 + 3.4 kcal/mole. 


It was assumed in these calculations that the heat 
liberated in reaction of unit mass of contaminants during 
formation of the alloy is the same as the heat liberated in 
the reaction of unit mass of the component substances. The 
error resulting from this assumption was estimated to be 
1 kcal/mole; this has been allowed for in determining the 
uncertainty in the heat of formation. 


The authors conclude by expressing their thanks to E, S. Smotritska for preparing the alloy, to A. N. Elistratova 
for her aid in carrying out the measurements, and to M. I. Ivanov for valuable advice. 
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MEASUREMENT OF THE ABSOLUTE CONCENTRATION 
OF THORON IN THE AIR IN INDUSTRIAL ENCLOSURES 


Yu. N. Burmistenko 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 505-506, 
December, 1960 
Original article submitted June 12, 1960 


Owing to its small half-life (54.5 sec), the concentration of thoron is usually estimated by an air flow method. 
The magnitude of the ionization current due to the thoron in the ionization chamber is a function of the speed 8 the 
air passing through the ionization chamber; and hence in order to determine the thoron concentration, the relation 
between the ionization current and the air speed (for constant concentration of thoron in the air) must be known. 


Let Qo be the concentration of Tn in air (disint /sec-1),V the volume of the ionization chamber (1 ), and w 
the speed of air (1/sec), The number of thoron atoms in one liter of air is then given by No = Q)/A where 2 is the de- 
cay constant of thoron. 


The number of thoron atoms entering the chamber per second is Now, while the number of thoron atoms leaving 
the chamber per second is Nowe 7) tt follows that the number of thoron atoms which disintegrate per second in 
the chamber is given by: 

a5 Qo set Narr 
N=N, (i—e ie (iA—e  ®). 


The ionization current is proportional to the number of disintegrations taking place in the chamber per second 
and the ionizing power of the a particles emitted during the disintegration of thoron. The current is given by: 


5 5 b (a) A - 
1=a-3,74-108e N=5,96-10-Ha Lo (1—-e. .® ata; (1) 
where a is a constant for the given chamber, 3.74 x 10° is the number of ion pairs produced by the a particles from 
Tn + ThA per disintegration, and e is the electronic charge. 


It is clear from this formula that as air speed increa- 
ses the ionization current tends to a limiting value which 
can be found with the aid of the l'Hospital rule for w> © 
and is given by 


lim Ty ,o0 = 5,96: 10-!4aQpV =kQ,V, (2) 
where k = 5.96 - 107 “a, 


The figure shows the ionization current as a function 
of the air speed for different volumes of the ionization 
chamber, The curves were calculated using Eq. 1 for lim 
Iya o= 1. 


The apparatus used to measure the thoron concentra- 
tion incorporates an ionization chamber, an electrometer, 
a flowmeter, and an airblower (vacuum main), 


The instrument is conveniently calibrated with the 
aid of a liquid radium standard. A correction was intro- 
duced for the difference in the ionizing power of a par- 
ticles from Rn and Tn (and their decay products), 


Ionization current due to thoron as a function of air 
speed, 
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The correction coefficient has been calculated by V. I. Baranov [1, 2] and is equal to 1.47. It follows that the 
conversion coefficient for the instrument is given by 


Qrn 
Rn 


Itn >= 


1,47 C/dw/min, (3) 


where Qpy is the amount of radon in the standard and Ipy is the standard ionization current. 


A liquid thorium standard cannot be used to calibrate the instrument,since the rate of production of the emana- 
tion by the standard would depend on the air speed and the construction of the bubbler. 


In order to determine the thoron concentration in air one must measure the ionization current for a constant 
speed of the air passing through the chamber. The thoron concentration can then be calculated from the formula 


Tent rn 1 : 
Se ee eG 
Qo V = /liter, (4) 


where Ip is the ionization current due to thoron, V is the volume of the ionization chamber, and 1/n is a correction 
which is determined from the curves given in the figure for a given air speed. 


The measured concentration of thoron may be affected by external y rays, or by the presence of radon in air. 
For this reason the ionization current due to thoron was determined as the difference between the ionization current 
measured while the air was passing through the chamber (Ig) and the ionization current measured with stationary air 
(1,), the difference being given by 


I rere bel ig (5) 
An exposure of 5 to 6 minutes is necessary for a virtually complete disintegration of the thoron. 
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ON THE FORM OF THE ABSORPTION CURVE FOR BETA-RADIATION 


IN ALUMINUM 


N. E. Tsvetaeva 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 507-508, 
December, 1960 
Original article submitted May 7, 1960 


The absorption curve for 8 radiation is often qualitatively described by the exponential oa I= heh 
where the absorption coefficient p is a function of the energy Eo of the B radiation. However, this representation of 
the absorption curve is very approximate.* Experimental data indicate that the form of this curve depends on ne 
geometrical conditions under which the measurements are carried out and differes appreciably from the exponential 
form [1, 2]. If the absorber is located next to the window of the counter, while the specimen is at a distance of sev- 
eral centimeters from it, then the absorption curve is very similar in form to the exponential. However, the initial 
part of the curve differs considerably from that predicted by the exponential formula. Some authors ascribe this to 
extraneous effects [2, 3]. This tendency to see an exponential curve where it does not really exist leads to consider- 
able errors in the determination of the absolute activity of 8 radiation with end-window counters, and also in the 
analysis of complex B radiation by the absorption method. 


Half-Layer Values of Absorption 


ters with atomic numbers Z = 4, 20,and 80. The computation was 

based on the data obtained by Seliger [4] for the absorption of mono- 

chromatic electrons, and the theoretical allowed B spectra given by 

Marshall in [5]. The computed curves correspond to the case where 

Cm) 70 et 80 ae 40 ae SO COLO the absorber is located next to the window of the counter and a narrow 
Thickness of Al absorber, mg/cm? beam of B particles is employed. 


J 
3 
us LSE Rr ea ee Ee eae d Data of 
f Rg Pe ed ee Ep, kev Zz 1/ai? 1/em? 
ASC ees ees es mg/em”_| mg/cm’ | [6] 
6 
5 SCE Le Eee a el 
; [a ee 4 2,80 4,5 if: 
(QUES a a a mo | mo | 2 | 45 | 
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2 i 
\ 4 19 26 
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1 22 aa ea coe ee ee 80 8 23 
pA 
eae Calera bares oe 
it. ES ee ee ee 
& ES Ge Ee ee We have computed by graphical integration the absorption curves 
3 [LONG se ee for simple B spectra with energies Ey = 200 and 600 kev, and Bemit- 
Biases | 


Bien 


Fig. 1. Computed absorption curves for 8 


Figure 1 shows that the computed absorption curves differ from 
radiation in Al, 


the exponential for small values of the absorber thickness. 


The table gives half-layer values obtained from the computed curves for the initial (di); and middle (d 1)m 
parts of the absorption curve, As can be seen from the table the values of (di); and (di) are very different. 
Moreover, (d1); depends on the atomic number Z of the B emitter. 


* The original absorption cur ves were based on ionization chamber data, These were closer to the exponential form 
than the geiger counter data. 
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Fig. 2. Computed (d4);= f(Ep) for 8 emitters with Z = 4, 20 
and 80. The points represent the values of (d1); obtained expe- 
rimentally [7]. 


Even in the case of the thin window, the use of (d1),, in correcting for the absorption of 8 radiation will lead 
to a large error in the absolute 8 activity. Thus, with a window thickness of 2mg/ cm’ the error is of the order of 
30%, In such cases one should use the function (41); = (Ep ) and take into account the atomic number of the B ‘emit- 
ter (Fig. 2). Figure 2 also shows some experimental points representing half-value layers for S$? Gat Go.sandsste, 
[7]. As can be seen these points agree with the calculations to within experimental error, except for s*°_ In the lat- 
ter case the low value of the energy of the 8 radiation (Ep = 0.167 Mev) made it difficult to determine (d1)j. 
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NEWS OF SCIENCE AND TECHNOLOGY 


SECOND ALL- UNION CONFERENCE ON NUCLEAR REACTIONS 
AT LOW AND MEDIUM ENERGIES 


A. A. Ogloblin and V. I. Chuev 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 509-529, 
December, 1960 


The Second All-Union Conference on Nuclear Reactions at Low and Medium Energies convened a othe 
July 21-28, 1960. Approximately 600 colleagues from nuclear research laboratories throughout the poet Union 
and scientists from 16 different countries were in attendance. The agenda of the Conference dealt with the most pres- 
sing problems in nuclear physics. 20 review papers were presented by Soviet scientists, outlining findings from over 
100 research projects, Reports by foreign delegates were also heard. 


The first day's sessions were devoted to research on few-nucleon systems, The lighter nuclei are of particular 
interest for research on nuclear phenomena. The study of the lightest nuclei may shed important light on such aspects 
of nuclear forces as are not manifested in nucleon-nucleon collisions (presence of multiple forces, etc.), Furthermore, 
data on nuclear reactions between the lightest nuclei are essential in the study of thermonuclear processes. 


A review of experimental research on few-nucleon systems was offered in a paper by N. A. Vlasov. Particularly 
detailed attention was given by the reported to an analysis of the spectra of neutrons and protons formed in reactions 
between H}, H’, H’, He’®, and He“ nuclei, These reactions are three- particle reactions, as a rule, and the character 
of the spectra of recorded particles depends strongly on interactions in the final state. Information on nuclear inter- 
actions produced in such experiments is not obtainable in many cases by scattering experiments. The repories also 
gave an account of experiments on measurement of neutron polarization in the T (p, n) He® and T (d>n)iHes reac- 
tions, and presented a vivid picture of the problem, now under heated discussion, of the possible existence of super- 
heavy isotopes of hydrogen and helium (H®, He®) and of the tetraneutron. 


Theoretical questions related to the structure of light nuclei were dealt with in the report by A. I. Baz". One 
of the chief defects in the shell model is the incorrect description of the periphery of the nucleus, failing to take 
into account formation of groupings of three and four nucleons, which have a high probability of appearing at the sur- 
face of the nucleus. The reporter dwelt in detailed fashion on the role of the nuclear surface in explaining the exis- 
tence of levels with large reduced widths. He also considered the influence of surface effects on the breakdown of 
charge invariance, observed in some individual cases, manifested in the nonconservation of isotopic spin in some reac- 
tions and in the appearance, in some nuclei, of levels bearing no analogy to the mirror nucleus (e.g., Li® and B®, Be® 
and B’, etc). 


H. Tyren (Sweden) told of experiments which he conducted in research on scattering of 450-Mev protons on D?, 
He’, Li’, and C” nuclei.. At this energy, quasi-elastic scattering of the incident protons on protons composing the tar- 
get nucleus takes place. An analysis of the experimental findings supports certain inferences as to the character of 
the momentum distribution of the protons in the nuclides studied. These findings are in harmony with shell-model 


predictions. I, S. Shapiro told of information obtainable on the structure of light nuclei through studying muon cap- 
ture by those nuclei. 


A host of data has been accumulated in recent years on the properties of nuclei by means of direct nuclear reac- 
tions going to completion without any compound nucleus forming in the interim. Studies of direct reactions, parti- 
cularly the stripping and pick-up reactions, have turned out to be highly valuable in nuclear spectrometry, since this 
research reveals information not only on spins and parities, but also on configurations of nuclear levels. Direct pro- 
cesses were discussed in review papers by A. A. Ogloblin, who surveyed the experimental work, and V. G. Neudachnik 
who discussed related theoretical problems. In the papers presented to the Conference, primary attention focused on 
a study of the structure of the nucleus in pick-up reactions, which at yet have received insufficient study. The (d, t) 
pick-up reaction on several nuclei ranging from Li to Al and several medium and heavy nuclei was treated in papers 
by N. A. Vlasov, S, P. Kalinin, A. A. Ogloblin, and V. I. Chuev. It was shown that excitation of hold levels is the 
predominant mechanism in pick-up reactions, leading to the determination of single-particle states inside the nucleus 
and determination of binding energy of the nucleon in those states. This sort of data was obtained for o%, F, and Zr 
nuclides. Mixing of the s- and d-states of the neutron in nuclei with A = 18 to A = 27 was also studied, as well as 
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the rotational structure of the Mg” nucleus. The other pick-up reaction (n, d) was studied by G. E. Velyukhov and 
A. N. Prokof'ev on Foe Ne”, po Se and by Laura Colli (Italy) and colleagues on nuclei from Se® to His, Findings 
on states of the proton in these nuclei were also obtained. 


In both review papers, emphasis was laid on the fact that subsequent application of the method of distorted waves 
with the use of optical potentials to stripping and pick-up reactions brought about remarkable agreement between ex- 
perimental data and theoretical predictions. For a correct choice of potential parameters, great importance is attach- 
ed to the data on polarization of particles as measured in the C!? (d, n) N® reaction by I. I. Levintov et al,, and by 
J. Nurzynskii and associates (Poland). E. Almqvist (Canada) told of studies on the mechanism involved in the Si® (d, 
p) Si® reaction with the aid of experiments in which the angular p-y correlations were measured and the variation 
in angular distributions of protons with energy was studied, 


Several papers were devoted to studies of the mechanism of (a, p), (p, &), (d, &), (@, t), etc. reactions. The 
nature of the angular distribution showed that these reactions proceed as a direct process of rather complex character. 
In particular, N. A. Vlasov and associates showed in their work that the (a, t) reaction is not equivalent to the (d, n) 
stripping reaction. The angular distributions in the (d, ~) reaction, as demonstrated in a paper by P. Cerineo (Yugos- 
lavia), evidence somewhat better agreement with the assumption that the reaction proceeds by a pick-up mechanism 
rather than by a knock-out process. 


An investigation of several direct processes occurring in response to bombardment of the c™ nucleus by protons 
at 150 Mev was discovered in a detailed paper submitted by P. Radvan (France). 


An investigation of scattering of nucleons on nuclei played a substantial role in creating the cloudy crystal ball 
model of the nucleus. A survey of research work on elastic scattering of protons was made by A. P, Klyucharev. This 
research encompassed study of proton scattering on targets of separated nuclides in the range extending from Z = 20 
to Z = 30 at proton energies below the potential barrier top. It was demonstrated that the angular distributions differ 
substantially for nuclei of even and odd A, if the numberof nucleons is close to a magic number, This is apparently 
related to the fact that an odd nucleon opens up the nuclear surface and appreciably increases the probability of ab- 
sorption of an impinging particle. With increased energy, the pattern of elastic scattering approximates to a diffrac- 
tion pattern and the isotope effect evokes a milder response from the shape of the angular distributions. A similar 
paper dealing with Cr, Co, Ni, Cu, and Zn nuclei was submitted by Chou Teh- Ming and associates (China). 


The present state of the cloudy crystal ball model of the nucleus was discussed in a paper by P. E. Nemirovskii, 
who reviewed work on inelastic scattering of nucleons. The reporter laid emphasis on the fact that the current inten- 
sive search for improved parameters for each separate nuclide is in contradiction to the fundamental concept of this 
model, Most of the papers also fail to provide any justification for the assertion that concentration of absorption at 
the surface of the nucleus leads to better agreement with experiment than a potential with an imaginary part distri- 
buted over the nuclear volume. N. Hintz (USA) adduced some empirical data at variance with this belief. Nemirov- 
skii devoted much attention to various applications of the cloudy crystal ball model to the theory of radiative capture, 
inelastic scattering, and polarization. 


J. Peterson (USA) gave an account of a new interpretation of giant resonances in total neutron cross sections, 
similar to the interpretation of the Ramsauer effect with reference to electrons, E, Almqvist gave a report on a study 
of giant resonances in total neutron sections, for nuclei from Mg to Zn. A resonance of this type was discovered for 
Ni®, 

A. Strzatkowski (Poland) delivered a report on elastic scattering of deuterons by some light nuclei. P. Heuge- 
louw (Netherlands) reported on the results of an investigation of elastic scattering of deuterons by Ag and Au. An ana- 
lysis of inelastic proton scattering data and (n, p), (, 2n), (n, pn) reactions was made, in his report, by U. Facchini 
(Italy). He pointed out that statistical theory is capable of explaining most of the efffects observed in such reactions, 
for the energy range 13-18 Mev. 


One of the sessions of the Conference was devoted to a study of nuclei with the aid of multiply charged ions, A 
review of work on Coulomb excitations of nuclei was made by I. Kh. Lemberg. He noted that the investigation of 
Coulombic excitation has been extended to cover light and medium nuclei, for which it is possible to obtain data on 
the lifetimes of excited states. It is interesting that data on the lifetime of Ne levels were obtained as a result of ex- 
citation of states in the bombarding nucleus, rather than in the target nucleus. 


1059 


A, Zucker gave a report on heavy-ion research in the USA. He noted that one of the most characteristic fea- 
tures of reactions with heavy ions at 10 Mev/nucleon energy is the important role of large angular Sea: Forma-~ 
tion of a compound nucleus is therefore not a basic process, and a neutron transfer reaction of the (N~, N~) type has 


a high probability of occurrence. 


E, Almqvist reported on heavy-ion experiments conducted in Canada. Scattering cross sections and reaction 
cross sections in the case of interaction between two carbon nuclei have disclosed some sharp resonances with an ano- 
malous high probability of decay of the compound nucleus Mg” into two C™ nuclei. An analysis of the findings leads 
to the conclusion that separate molecular states exist in Mg” with lifetimes of 107” to 10° sec, which may be rep- 
resented as two oscillating c” nuclei. 


A problem of enormous interest in modern physics is the effect of resonant absorption of gamma rays without 
recoil,dicovered by R. Mossbauer in 1958, The essence of the Mossbauer effect is that the recoil momentum of a 
gamma photon emitted by a nucleus under certain conditions is imparted to the surrounding crystal as a whole, rather 
than to the nucleus alone. The gamma photon thereby carries off all the energy of decay and may be resonantly ab- 
sorbed by a nucleus similar to the original one. This phenomenon and its numerous applications were discussed in a 
report by R. Mossbauer (West Germany), heard with intense interest. This report laid the foundations of the theory 
underlying this effect, and described experiments in which the effect was observed. Extremely narrow line widths (as 


narrow as 10 °° of the nuclear transition energy for Zn°’) have been brought within range of measurement by the Moss- 
bauer effect. 


The Mossbauer effect opens up extensive opportunities for studying the nature of intraatomic magnetic and elec- 
tric fields, the magnetic and electric moments of excited states of nuclei, the effect of the detailed structure of elec- 
tronic shells on the position of nuclear levels, R. Pound (USA) reported on a novel utilization of the Mossbauer effect 
for testing the general theory of relativity. He observed a change in the frequency of photons under terrestrial condi- 
tions as the photons progressed in a gravitational field. The measured frequency difference of upward and downward 
moving photons agrees with predictions of relativity theory to within 4%, 


A report by F. L. Shapiro discussed the classical interpretation of the Mossbauer effect. This reporter also sur- 
veyed experimental research performed in the USSR on the study of this interesting phenomenon. 


In the area of photonuclear reactions, primary attention was devoted to research on the shape of the giant re- 
sonance curve. A review of experimental work in this direction was made by L. E. Lazareva. The study of the reso- 
nance curve was made for oxygen. Measurements in the 21-23 Mev range of energies were carried out with high re- 
solution, and the position of the narrow maximum detected in absorption cross sections was found to be in excellent 
agreement with theory. Some interesting data on the mechanism responsible for photodisintegration of nuclei were 
reported in papers by R. M. Osokina and associates. They showed that protons and neutrons are emitted as a result of 
direct interaction. In a paper by O. Bogdankevich and associates, the inelastic scattering cross section for gamma 
photons was measured and some features near the threshold of the (y, n) reaction were revealed. 


A detailed report on the shape of the giant resonance curve for nuclei with large deformations was delivered 
by E. Fuller (USA). The giant resonance peak splits into two humps in the case of deformed nuclei. The deforma- 
tions and quadrupole moments of the nuclei may be found from the energy relations of both resonance peaks, and 
turn out to concur closely with the values arrived at in other experiments. Data on elastic scattering of gamma 
quanta by Ta, Ho, and Er indicate the existence of nuclear combination scattering. 


A review of theoretical research on the mechanism responsible for photonuclear reactions was presented by 
A. M. Baldin. 


The general regularities evident in the spectra of gamma photons produced in radiative capture of slow neutrons 
were discussed in a paper submitted by L. V. Groshev. At the present time, spectra for most elements have been de- 
tived, and it has been found that the shape of the spectrum is in large measure determined by the distribution of level 
densities. The development of new gamma-ray spectrometers with high resolving powers has in some cases enabled 
observation of excitation of single-particle and collective states. This was illustrated by the example of the nuclide 
ca™, whose spectrum was measured in the paper presented by L. V. Groshev and associates, 


S. P, Tsytko presented a paper on radiative capture of protons. 


Models of the nucleus were discussed at the last two sessions. A report on the shell model was presented by 
V. V. Balashov. The reporter emphasized that the shell model is in its present state (where residual interactions, 
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etc, are taken into account) capable of accounting for a host of experimental facts, particularly in the field of light 
nuclei, 

A. S, Davydov laid great stress, in his report, on the model of a deformed nucleus lacking axial symmetry, a 
concept developed by himself and co-workers. The presently available empirical data on positions of nuclear levels, 
moments, and transition probabilities are not at variance with this model. D. P. Grechukhin told of one possibility 
for making a choice between the axial and nonaxial symmetry models of the nucleus. 


Yu. T. Grin" presented a paper devoted to the application of the recently developed superconductivity theory 
to the nucleus. The reporter gave an account of the theoretical work carried out in the recent period by A. B. Migdal, 
S. T. Belyaev, S. I. Drozdov, D. F, Zaretskii, Yu. T. Grin', and others in which the effect of paired correlation on 
the various properties of the nucleus comes to the fore, L. A. Sliv took the floor on the same problem. 


F, Villars (USA) dwelt upon several problems related to the theory of strongly deformed nuclei, particularly on 
quadrupole oscillations. G. Brown (United Kingdom) told of the results of calculations of nuclear-shell wave func- 
tions, with interactions between a nucleon and a hole taken into account. 


This Conference on nuclear reactions undertook a much wider range and scope of problems than its name would 
indicate. The addition of several new problems to the Conference agenda made it possible to outline the basic trends 
evident in the development of low-energy and medium-energy nuclear physics, and to shed light on the most inte- 
resting results arrived at in recent work. 


The materials of the Conference will be published by the Academy of Sciences of the USSR as a separate sym- 
posium, during the second quarter of 1961. 
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INTERNATIONAL CONFERENCE ON PLUTONIUM METALLURGY 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 509-529, 
December, 1960 


The International Conference on Plutonium Metallurgy, organized by the Commissariat de 1'Energie Atomique 
and the French Metallurgical Society, convened at Grenoble (France), April 19 to 22, 1960. 200 persons from 20 
countries were in attendance at the conference, to hear over 40 papers submitted prior to the conference. 


The first session was devoted to the physical properties and the metallurgy of plutonium. It may be said that 
six allotropic modifications of plutonium and their properties are known at the present time (properties such as elec- 
trical conductivity, emf, specific heat, heat of transformation, magnetic behavior) to some significant degree. The 
anomaly of electrical resistance at 100°K is yet to be explained, the crystalline structure of the a and B phases is 
still obscure, and the problem of the negative coefficient of expansion of the 6 phase has not been fully resolved. 


In a report by Kerjean (France), a description was given of the technology of fabrication of 200 g plutonium 
ingots via reduction of calcium tetrafluoride. 40% excess of calcium is taken (the British prefer a somewhat smaller 
excess); the burden is compacted into pellets. Twelve pellets are charged semiautomatically into a calcium fluoride 
crucible, which is high-frequency heated in an argon atmosphere under the conditions; heating to 600-800°C for 1 
min 30 sec, leaving to stand for 45 sec, 97.3% metal yield. All of the steps are carried out in a line of argon-filled 
glove boxes. Since plutonium tetrafluoride is an intense source of fast neutrons, extreme precautions must be observed 
when handling it (keeping a distance of 80 cm from the source, avoiding handling of the substance more than 2 to 3 
hours daily). 


A variety of methods for producing pure plutonium received some measure of attention, in particular electro- 
lysis with careful prior purification of the original plutonium nitrate. Calcium thermal reduction in crucibles of 
magnesium oxide resulted in 99.8% pure metal at Los Alamos, 


The second session was devoted to plutonium alloys. Phase diagrams of the Pu — Nb, Pu—Cu, Pu —C systems 
were reported on, and results were presented for research on ternary systems, especially the Pu — U — Mo system, 
which has received detailed study with respect to the kinetics of phase transformations. 


The third session heard reports of findings on the ceramics technology of oxides of plutonium and PuO,-UO, 
systems, in connection with the great interest manifested in the use of these materials as nuclear fuel. Several va- 
riants of the sintering process for PuO, were examined, and it was shown that sintering of PuO, follows a path similar 
to that for UO2. Work on sintering studies, behavior under irradiation, and results of measurements of the thermal 
expansion coefficient for PuO, and PuO2-UO, solid solutions were reported on by representatives of the USA. 


At the concluding session, problems involving applications for plutonium in thermal reactors (fuel elements in- 
corporating Pu — Al alloys) and in fast breeder reactors (the LAMPRE reactor with a liquid eutectic for Pu — Fe melt- 
ing at 410°C, and the Rhapsodie reactor with fuel elements of Pu — U — Mo alloy) were discussed. The mechanical 
properties of plutonium and problems involving the compatibility of plutonium alloys and cladding materials (stain- 
less steel, niobium, etc.) were examined. 


The results of current investigations of plutonium-base nuclear fuel may be summarized as follows. 


Pure plutonium is unsuitable as nuclear fuel because of the presence of allotropic transformations and other un- 
favorable properties. 


Liquid plutonium-base fuels are quite expensive. 


Plutonium carbides have not been investigated to any great extent, and it is still not proved that their thermal 
conductivities and irradiation behavior are the same as those of uranium carbide. 


Ternary alloys have excellent properties, but are difficult to clad. 


The most favorable results have been obtained with plutonium oxides and PuO, — UO, solid solutions, despite 
some unresolved problems which are still plaguing the researchers, 


1062 


Ae DU ESLATVIANSINSTITUTEVOh PHYSICS 
Yu. Koryakin 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp, 512-514, 
December, 1960 


At Salaspils, 20 km from Riga, construction work is being completed on the first nuclear research reactor in 
the Baltic republics, the IRT reactor of the Institute of Physics of the Academy of Sciences of the Latvian SSR (Fig. 1). 
This reactor, IRT-type, is rated at 2000 kw(th). Similar reactors are being built at various points throughout the 
Soviet Union. 


At the present time, finish-up work and assembling is going on in the reactor hall (Fig. 2). Some of the admi- 
nistrative quarters are already completed, and work has been started in them on immediate preparation of research 
work to be carried out on the reactor. The control and measuring instruments, reactor control system, shielding, and 
dosimetry laboratory will soon be operational. The reactor tank has been mounted in place and the heavy concrete 
biological shielding is nearing completion. The reactor assembling job and placing of instrumentations were carried 
out under the supervision and with the direct participation of workers attached to the Institute's reactor laboratory. 
Colleagues of the I. V. Kurchatov Order of Lenin Institute of Atomic Energy and other institutions were also of great 
assistance to them in this tank. 


Many of the contracts and operations involved in erection of the reactor were carried out through the efforts of 
the Institute of Physics and enterprises under the jurisdiction of the Latvian Council of the National Economy. For 
example, the reactor control desk was designed under the supervision of colleagues of the reactor laboratory of the 
Institute, and was fabricated in Latvian plants. 


Builders performed a number of construction tasks which were new to their experience, in the course of erecting 
the reactor plant. For example, a method for preparing and pouring heavy concrete was put into practice. In the 
zone of highest radiation, heavy concrete of 6.5 g/ cm? density was used, while concrete of 5.2 g/ cm? density was 
used in other shielding sections. 


In line with research schedules projected in advance, some modifications were introduced into the standard 
reference design. For example, the single hot cell originally envisaged was extended to two hot cells before plans 
were frozen, 


The progress of work on the reactor project was a constant concern to the government of the Latvian SSR, par- 
ticularly to J. V. Peive, Chairman of the Council of Ministers, who is a frequent visitor to the reactor site. 


The start-up of the reactor will mean a significant boost to the activities of the Institute of Physics in nuclear 
physics problems, and will encourage penetration of the achievements of nuclear energy uses into the national eco- 
nomy of the Baltic republics. Scientists from Lithuania and Estonia will also participate in research work associated 


with the reactor. 


The following trends were outlined for the work of the Institute at the August 1960 session of the Presidium of 
the Academy of Sciences of the Latvian SSR: physics of the nucleus and uses of isotopes, solid state physics, and 


magnetic phenomena. 
The avenues of research involved in the first category are, with particular emphasis, the following: 


1) nuclear spectroscopy of short-lived nuclides and capture gammas. Properties of radioactive nuclides (emis- 
sion spectrum, lifetime, etc.) with half-lives in the second and minute range will be investigated; 


2) investigation of radiations by means of magnetic spectrometers. Methods will be developed for research on 
circular polarization of gamma quanta in response to capture of polarized neutrons, and the quantum characteristics 
of the energy levels of compound nuclei will be determined. Magnetic spectrometers will be relied upon for studying 
the spectra of internal-conversion electrons of radioactive nuclides with half-lives ranging from several seconds up to 
two to five days. An electromagnetic dispatching system will be built for delivery of irradiated specimens from the 
reactor to the spectrometer for similar research on radioactive nuclides with half-lives in the second and minute range. 
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3) studies of neutron fluxes in the core and experimental channels of the reactor; 


4) theoretical research on ionization and excitation of atoms and ions by slow electrons, This research is dic- 
tated by the needs of applied spectroscopy, astrophysics, and other branches of physics and chemistry with regard to 
derivation of the fundamental characteristics of atoms and molecules, particularly and in the first instance for data on 
the effective excitation and ionization cross sections of atoms and molecules, In particular, a-research project is in 
the hopper for computing the excitation functions of various atoms by programming them into an electronic computer; 


5) the study of transient response in relay-type radioactive sensors. These investigations are called for by the 
fact that the fundamental parameters of relay-action instruments (sensitivity, speed of response, reliability, perfor- 
mance stability) are sensibly dependent on the radiation detectors used and on the circuitry of the relay recorders. 
The interrelatedness of these parameters is felt particularly where transients are involved; the transients are determined 
by the changes in radiation intensity in the various recording processes. 


Fig. 1. General view of the reactor building under construction at the Salaspils reactor site. 


Another of the laboratories attached to the Institute of Physics, one engaged in studies of radioative automation 
techniques, has already started work on its program for developing techniques for utilizing radioactive isotopes in 
science and technology. This program includes a study of the effect of random statistical processes on basic instru- 
ment constants in radioactivity measuring instruments, investigation of the yield and energy spectrum of bremsstrah- 
lung and characteristic radiation from beta sources, evolving efficient recording techniques, developing new sources 


of a SO rays and new beta preparations, a scientific development of fundamentals of standardization and nor- 
malization of process control instrumentation using nuclear radiations 


The topics included in the second category of research mentioned above are: 


1) the study of problems of radiation effects on the magnetic and electrical properties of metals and alloys; 


2) investigation of nuclear radiation effects on semiconductors and dielectrics; 
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3) semiconductor diffusion studies; 
4) nuclear magnetic resonance studies of the structure of matter; 


5) investigation of magnetohydrodynamic phenomena in flows of liquid metals at low Reynolds numbers; 


6) study of techniques of electromagnetic transport of liquid metals. 


The Latvian Institute of Physics has been engaged for the last five years in research on developing the founda- 
tions of electromagnetic pump theory. During that time, an approximate theory and calculation procedure has been 
arrived at applying to the basic variants of electromagnetic pumps, induction-type and conduction-type pumps. Sev- 
eral experimental pump designs have been worked out, with excellent performance characteristics. In 1960, work was 
started under contract with the Leningrad Council of the National Economy on building electromagnetic metering 


pumps for feeding pressure casting machines, The savings realizable from putting such metering pump units into pro- 
duction practice run into tens of millions of rubles annually. 


It would be no exaggeration to state that a whole 
school of specialists on electromagnetic pumps, a line of 
equipment with widespread applications in nuclear industry 
and research, has been trained under the Latvian Institute 
of Physics. 


The scope of the problems facing the Institute is thus 
a true reflection of the enormous gains registered by Latvian 
physicists and scientists in the field of peaceful uses of nuc- 
lear energy. 


The building of the Salaspils research reactor marks 
only the first stage in the expansion of nuclear energy ap- 
plications in Latvia. The Salaspils reactor will provide the 
basis around which a nuclear research center is projected 
for the coming years. 


A building for physics and engineering work on radio- 
active materials and radiations and advanced techniques 
for use in the national economy is one goal. Another is the 
construction of a radiochemical building for the production, 
processing, and shaping of pile-produced radioisotopes. 


A building with cryogenic equipment for the produc- 
tion of liquid helium and liquid hydrogen will be built on 
the premises of the coming research centers, as well as a 
radiation circuit with high gamma radiation level required 
for work in radiological physics, radiochemistry, and radio- 
biology. 


Plans now being drawn up for this center include a 
Fig. 2. Laying the bottom of the reactor pressure vessel _ building housing a linear accelerator and another for a 25 
in place. Mev proton accelerator. It will be used to investigate be- 

tween charged particles and nuclei, and to produce some 
radioactive isotopes which are difficult to obtain from reactors. Finally, living quarters for the research workers at- 
tached to the center will be set up in the environs of the center. 


After completion of a whole complex of nuclear research buildings at Salaspils, the Latvian Institute of Physics 
will place the most advanced equipment at the disposition of the center, which will be capable of resolving problems 
of first-ranking interest to the national economy in the area of peaceful uses of atomic energy. 
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THE HUNGARIAN INDUSTRIAL EXHIBIT AT MOSCOW 


Yu. Mityaev 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 515-516, 
December, 1960 


The Gorky Central Park of Culture and Recreations was host during August and September, 1960 to the Hungarian 
Industrial Exhibit. Instruments reflecting the most highly developed branch on nuclear industry in Hungary, that of 
instrument design, figured prominently among the exhibits. The exhibits attested to the great gains achieved by Hungary 
in peaceful uses of atomic energy (Fig. 1), Instruments fabricated in the Republic several years back were placed on 
exhibit for comparison purposes, illustrating in sharp contrast what great strides forward Hungary has made in the field 
of nuclear instrument design. 


A clinical device of original design, a gamma-ray 
scintillation isotope scanner (Fig. 2) designed to detect and 
map cancerous and other tumors with the aid of radioactive 
isotopes, appeared on exhibit. A harmless dose of a radio- 
active preparation is introduced into the organism of the 
patient, in a form such that it will be deposited in the site 
of the supposed tumor (e.g., a preparation containing ee 
would be administered to locate a thyroid tumor). Then 
the isotope scanner automatically maps the radioactivity 
of the patient's body. Radiation maps taken with two dif- 
ferent body orientations reveal not only the location but 
the size as well of the tumor, The gamma-ray isotope scan- 
ner operates with a scintillation-counter head, and is com- 
pletely automatic. Different modes of performance may 
be resorted to in order to obtain low-definition or high-de- 
finition radioactivity maps. About 80 such instruments are 
available at present. All of them are being employed suc- 
cessfully in clinics and hospitals throughout the Hungarian 
Peoples Republic. 


A variety of instruments designed to measure radia- 
tion intensities under laboratory, plant, and field conditions 
alike were given prominent display at the exhibit. These 
are portable transistorized devices, with low weight and 
little bulk; they can be employed to great advantage in 
agriculture, biology, and uranium prospecting. Most of the 
instruments are designed to work from either a scintillation counter of a halide gas-discharge counter, Some of the 
portable instruments have a small built-in Co™ source for determining radiation levels by the comparison method. 


Fig. 1, At the Hungarian Industrial Exhibit. Photo by 
S, Blinov. 


Reliable pocket dosimeters (together with charging device) designed for dosimetric monitoring of personnel work- 
ing under exposure to x-rays and gamma rays were also in display. Among the other dosimetric instruments, there 
was an acoustic threshold indicator for signaling an alarm whenever a predetermined radiation level is surpassed. 


A set of instruments for precision measurement of radioactive speed specimens was put on display. The set con- 
sists of an original-design lead castle (Fig. 3) housing the test preparation, a broadband 70 Mc amplifier for a pulse- 


height analyzer with resolving time of 5 sec, and a power supplies unit. The design of the lead castl facilitates the 
study of both solid and liquid specimens. 


In addition, a varied display of radio engineering equipment was available: Scalers with wide bandwidht, high- 
voltage stabilized power supplied, timers and clocks for remote control of various physics instruments, standard mo- 
dular units for synthesizing measurements circuits, etc. Scintillation crystal of vinyltoluene produced in Hungary can 
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Fig. 2, Gamma-ray scintillation isotope 
scanner. Fig. 3. General-purpose lead column. 


be fabricated with ease into wafers of any shape, forming a substrate for deposition of a layer of some other scintil- 
lator, e.g., Nal. 


A characteristic feature of the exposition was the fact that the overwhelming majority of the instruments ex- 
hibited are being mass-produced by Hungarian industry and are being used in many applications not only in scientific 
research, but also in the various branches of industry, medicine, and agriculture. 
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AT THE JAPANESE INDUSTRIAL EXHIBIT IN MOSCOW 
Yu. Koryakin 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 516-517, 
December, 1960 


The Sokol'niki Park of Culture and Recreation was host, during August and September, 1960, to the Japanese 
Industrial Exhibit. On display was a broad variety of electronic equipment produced by Japanese firms. Of partic- 
ular interest were several types of electron microscopes of various magnifications (Fig. 1) and a large number of 
electronic devices by Yokogawa Ltd. Yokogawa instruments are primarily for process control and enjoy a broad 
range of application. 


sas Ay - a 
Fig, 2. Yokogaw [a complete electronic control system]. 


Fig. 1. Japanese manufactured electron microscope. 


An arrangement of compressed size for quanti- 
tative and qualitative x-ray analysis of samples and 
the composition of test materials attracted some 


attention (see Fig, 3), The facility is designed for 
metallurgical, chemical, textile, ceramic, and other Fig. 3, Setup for x-ray analysis of specimens and composi- 
branches of industry, tion of materials, 
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NEW REGULATIONS GOVERNING HANDLING OF RADIOACTIVE 
MATERIALS AND SOURCES OF IONIZING RADIATIONS 


L, Lazareva, N. Leshchinskii, P. Moiseitsev, 
V. Sinitsyn, and A, Shtan' 


Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 525-529, 
December, 1960 


The State Committee on the Uses of Atomic Energy under the jurisdiction of the Council of Ministers of the 
USSR and the State Health Inspection Commission of the Ministry of Public Health of the USSR have evolved, and 
promulgated on June 23, 1960, a new set of “Health and Safety Regulations Governing Handling of Radioactive Mate- 
rials and Sources of Ionizing Radiations," No. 333-60 [1], which now replace the previous "Health and Safety Regula- 
tions Governing the Transportation, Storage, and Handling of Radioactive Materials," No, 233-57 [2] promulgated in 
1957, 


The regulations will be communicated to all institutions, enterprises, laboratories, and other organizations un- 
der all ministries and local authorities engaged in transporting and storing radioactive materials and sources of ioniz- 
ing radiations, They are obligatory regulations governing the design, construction, reequipping, and operation of en- 
terprises, institutions, laboratories, and other buildings and facilities intended for work with radioactive materials and 
sources of ionizing radiations. 


The new regulations take due account of experience gained in using radioactive materials and ionizing radia- 
tions in the Soviet Union during the recent period, as well as recent data on their effects of living organisms and on 
materials, and the basic recommendations of international bodies relevant to these questions [3, 4]. 


New critical tolerance levels for ionizing radiations and new critical tolerance concentrations for radioactive 
substances in water and air were put into effect at the same time as the regulations referred to. Other new standards 
include a unique radiation hazard symbol for proper marking of objects, rooms, equipment, instruments, materials, 
etc., inside or outside of which a radiation hazard is present, nameplates for enterprises, institutions, laboratories, etc. 
authorized to handle radioactive materials, and several other regulations, requirements, statutes and forms. 


The regulations are based on the need to ensure safety in handling of radioactive materials, and incorporate 
the latest critical tolerance dosage. These regulations stipulate 5 rem annually as the permissible dosage for any wor- 
kers exposed to radioactive materials or sources of ionizing radiations, with reference to all forms of external irradia- 
tion. This form corresponds to 0.1 rem weekly for personnel continually exposed, For personnel working with a radius 
of exposure to the radiations, but not directly engaged in handling radioactive materials or sources, the yearly dosage 


is set at 0.5 rem. 


The regulations are careful to strictly delineate requirements governing handling of radioactive materials in 
properly closed containers, i.e., in such a state or packed such that the radiations do not propagate to the surrounding 
medium, so that they can act only as external radiation sources, distinguishing them from radioactive substances pre- 
sent in open form whereby the radiations are not hindered from propagating into the surrounding medium and affect- 
ing the organism both from outside and by penetration into the organism. To provide safety in handling of radioac- 
tive materials present in closed form and in handling and using other sources of ionizing radiations (x-ray tubes, ac- 
celerators, etc.), only biological shielding need be specified; i.e., attenuation of the radiations is required at those 
points where personnel may suffer exposure, to keep the dosage well below critical levels. 


The regulations spell out concrete measures to be fulfilled in building and using high-level isotope facilities, 
x-ray diagnostic equipment, radiotherapeutic equipment, nondestructive radiography facilities, and other equipment 
and instruments involving the use of radioactive substances. In observance of these requirements, instruments may 
be employed without special restrictions in workshops, installations, laboratories, and other work places, thereby sig- 
nificantly broadening the field of possible applications without incurring a hazard of radiation injury. 
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Requirements and recommendations governing work with radioactive materials present in open ‘oo are spel 
out painstakingly in the regulations, All radioactive materials are classified according to their radiation eet) ae 
diation effects when radiation enters the organism) into four groups: Group A includes radioactive materials of Pate 
cularly high radiation toxicity (sree Pow’ Rae etc.), group B includes radioactive Faas pe nigh radiotoxicity 
level (Ca**, Co®™, 181, etc.), group C includes materials of average radiotoxicity (Na™, P**, Fe”, I", ete.), and group 
D embraces low-toxicity isotopes cH, c¥, cr®°4, etc.). Assigned to group A are, primarily, those radi qarrve substan- 
ces whosé critical concentrations in the air of rooms occupied by personnel is not more than 1 - 10” curie/liter. 
Group B has assigned to it substances whose critical concentrations range from 1 - 10°" to 1-107 curie/ liter, and 
group D includes substances for which concentrations greater than 1 - 107° curie/ liter in air are permitted. 


Depending on the degree of radiotoxicity and the activity of the radioactive substances at the working site (start- 
ing at 0.1 microcurie for the most highly radiotoxic substances), three classes of operations involving exposure to ra- 
diation have been singled out. Work operations assigned to class III may involve the use of ordinary, nonshielded work- 
ing places equipped in compliance with the regulations applying to chemical laboratories, and in some cases work 
with nonradioactive substances may be authorized for the same places. Work operations assigned to class II must be 
carried out in rooms specially set apart and,furthermore equipped with specialized or adapted equipment, depending 
on the nature of the work. Work operations assigned to class I must be restricted to separate work places with special 
zonal planning, safety partitions, specialized ducting and hoods, specialized equipment for high-level work (glove 
boxes, junior caves, etc). 


Institutions, enterprises, and laboratories are grouped into three categories depending on the annual needs for 
radioactive materials in open form, since the amount of hot wastes and contamination of air and drainage passages is 
largely dependent on the magnitude of those needs, 


The pertinent class and category are governing factors for institutions, enterprises, and laboratories with respect 
to specifications for planning, isolation, and equipment of working places and rooms, the system specified for waste 
collection and disposal, drainage, ventilation, the use of various types of operating equipment, etc. 


The regulations devote considerable space to problems concerning the correct organization of the production, 
transportation, inventorying, issuing, and storage of radioactive materials, since proper management is a deciding fac- 
tor in the preservation, safe handling, and rational use of such substances. 


The regulations forbid introduction of radioactive materials into foodstuffs. Applications for radioactive subs- 
tances in the various branches of the national economy via introduction of those substances into the finished industrial 
product is permissible under conditions that the residual activity (in curies/kg) of the end product not exceed by more 
than 100 times the critical tolerance concentration (in curie/ liter) for water in open reservoirs in the case of radio- 
active substances having half-lives up to 60 days. For substances with half-lives longer than 60 days, residual activity 
must not exceed 10 times the critical level for water in open reservoirs. 


Requirements governing the collection, disposal, and burial of radioactive wastes, including regulations cover- 
ing conveyances and burial sites, are listed in special subsections. 


Particular attention is devoted in the regulations to measures for individual shielding, personal hygiene, and do- 
simetric monitoring. A list of recommended protective clothing, "Don't Do. . .” rules for working with radioactive 


substances and sources of ionizing radiations, techniques and means for decontamination, and other useful information 
are appended to the regulations. 


An appendix of teminology and technical definitions is included to facilitate understanding of the regulations. 
Also appended are formulas and tables, an ample selection of examples illustrating radiation shielding calculations, 
of concentrations of hot substances, and other valuable information. 


The introduction of these new regulations remove any existing ambiguities in realm of applications of radio- 
active substances and sources of ionizing radiations, and establish the requirements and specifications corresponding 
to the present advanced state of knowledge about the effects of radioactive substances and ionizing radiations on the 
organism, which will contribute to the safe use of these tools on a broader scale. 
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Sanitarnye pravila raboty s radioaktivnymi veshchestvami i istochnikami ioniziruyushchikh izluchenii [Health 
and safety regulations governing the handling of radioactive materials and sources of ionizing radiations]. 
Moscow, State Atom Press (Gosatomizdat, 1960). 

Sanitarnye pravila perevozki, khraneniya, ucheta i raboty s radioaktivnymi veshchestvami [Health and safety 
regulations governing the conveying, storage, inventorization, and handling of radioactive substances]. Moscow, 
Maritime Industry (Minsudprom, 1957). 

Bezopasnoe obrashchenie s radioioztopami [Safe handling of hot isotopes]. Vienna, Internat'l. Atomic Energy 
Agency, 1988. 

Recommendations of the International Commission on Radiological Shielding. Munich, ICRS (1959). 
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NEW LITERATURE 


Books and Symposia 

Accelerators, A symposium. Moscow, Gosatomizdat, 1960 
124 pages, 3 rubles, 60 kopeks 
‘Translated from Atomnaya Energiya, Vol. 9, No. 6, pp. 530-531, 
December, 1960 


This symposium is composed of original articles decoted to advanced linear and cyclic accelerators. The first 
three articles discuss linear electron accelerators: A 6 Mev linear with constant phase velocity, installined at the Uk- 
rainian Physics and Engineering Institute, is described; some problems in particle dynamics are studied and the pro- 
cess of electron bunching in a linear electron accelerator is discussed. A new plan for particle extraction from a syn- 
chrocyclotron is proposed in the fourth chapter. The fifth chapter presents a technique for beam extraction from a 
betatron or synchrotron doughnut by unilateral displacement of the equilibrium electron orbit outside the operating 
region, a technique developed at Tomsk Polytechnic Institute. The sixth chapter describes a ferrite frequency trans- 
lation device for shifting a cyclotron into the synchrocyclotron mode, The seventh chapter takes up such problems as 

. shaping axisymmetric magnetic field configurations by positioning of annular shims, The eighthchapter deals with the 
generation of multiply charged ions in a cyclotron. The ninth chapter describes a periodic-field cyclotron for accel- 
erating multiply charged ions. The last three chapters discuss the effect of multiple scattering and emission in beam 
stacking, a technique for measuring the duration of an ion cluster in a cyclotron, and the results of an experiment on 
annular electron beam scanning at 3 Mev. 


N. D, Fedorov. The cyclotron. Moscow, Atomizdat, 1960 
88 pages, 2 rubles, 30 kopeks 


The brochure tells of the cyclotron's background, delving into the early days of the cyclic resonance in accel- 
erators, The operating principle underlying the cyclotron is discussed, and the problem of enhancing ion energy and 
beam intensity is gone into. Systems for extracting ions from the accelerator chamber are described. Basic informa- 
tion is given on cyclotron design, and a description of a cyclotron with an azimuthally varying magnetic field appears. 
The major applications of cyclotrons are indicated. 


The brochure is written for a broad public. 


M. A, Bak and Yu. F, Romanov. The neutron. Moscow, Atomizdat, 1960 
82 pages, 2 rubles 15 kopeks 


This brochure describes the discovery of the neutron and explains its fundamental properties. An account is 
given of interactions between neutrons and matter, fission of heavy nuclei in response to neutron bombardment, neu- 


tron sources, production of monoenergetic neutron beams, discovery of the antineutron. Neutron dosimetry and neu- 
tron shielding techniques are touched upon. 


D. N, Trifonov. Element 61, Its Past, Present, and Future, Moscow, Atomizdat, 1960, 
56 pages, 84 kopeks. 


This booklet is devoted to the element promethium, of atomic number 61. A popular-style account is given 
of the history of its discovery and methods of winning and purifying it, closely related to the history of the other rare- 
earth elements, one of the most interesting problems in modern inorganic chemistry. 


The booklet is aimed at the general reading public. 
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G. F, Silina, Yu, I, Zarembo, and L. E. Bertina, Chemical Technology and 
Metallurgy of Beryllium, Edited by V. I. Spitsyn, Moscow, Atomizdat, 1960, 
120 pages, 3 rubles, 60 kopeks. 


The basic problems involved in beryllium production and applications are discussed in this book, Chemical and 
physical properties of beryllium metal and beryllium compounds of technological interest are described. In the chap- 
ter devoted to industrial techniques for producing beryllium compounds, methods are given for detecting beryllium 
and producing technical-grade beryllia, purifying beryllium hydroxide to nuclear-grade specifications, producing am- 
monium fluoroberyllate, beryllium chloride, and beryllium carbide. Of various techniques developed to produce be- 


ryllium metal, those having some industrial significance, such as thermal reduction of halides, and fused-salts elec- 
trolysis, are discussed. 


Data published in the literature during the past 15 years are used in the book. 


The book is written for specialists working in the field of beryllium production and applications, and for students 
in metallurgical schools. 


V. V. Fomin, The Chemistry of Extractive Processes, 
Moscow, Gosatomizdat, 1960,168 pages, 6 rubles, 60 kopeks. 


This book, on the basis of published findings, summarizes current concepts of the extraction mechanism. The 
dependence of distribution constant on extraction conditions is discussed, and methods for evaluating this dependence 
are compared, Detailed numerical calculations are absent from the book; some mathematical derivations are left to 
the appendix, in order to simplify the presentation. A reference bibliography names 236 titles. 


The book is aimed at scientific workers and students interested in extraction problems. 


Yu. V. Sivintsev. Background Radiation of the Human Organism. Moscow, 
Gosatomizdat, 1960. 96 pages, 3 rubles. 


On the basis of experimental aata published in the literature, the components of external and internal irradia- 
tion of the human organism are discussed and the dose rates of the background irradiation of different parts of the hu- 
man body, such as the gonads, lungs, and skeleton are evaluated. A comparison is drawn between these doses and 
dangerous mutagenic and carcinogenic doses. Proof is presented that the background radiation dose may be accepted 
as a basis in judging maximum radiation tolerance levels. A list of the pertinent literature consulted appears at the 
end of the book, with over 200 titles mentioned. 


The book is written for physicists and biologists working in the field of research on biological effects of physical 
radiations, and for other persons interested in this problem. 
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